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ABSTRACT 


The Navy of tomorrow will require a robust and 
reconfigurable power system capable of supplying power not 
only to large high power propulsion loads but to growing 
combat system loads like high power radar and pulse loads 
such as rail guns and f ree-electron lasers. A critical 
component in such a system is the phase controlled 
rectifier. As such, the issues associated with the 
inclusion of a power electronics rectifier need to be 
addressed. These issues include input Alternating Current 
(AC) interface requirements, the output Direct Current (DC) 
load profile, and overall stability in the presence of non¬ 
linear loads. Understanding these issues and determining 
the means of assuring compatibility with a Navy all¬ 
electric ship is the focus of this thesis. 

By using a Simulink® model of a variable parameter 
load, several multiple-pulse count, high power rectifiers 
were exercised. The Simulink® results were compared to the 
linearized small signal transfer function analysis results. 

These experiments led to the conclusion that 
increasing the pulse count and output filtering reduces the 
input interface current distortion. However, there are 
tradeoffs in terms of complexity and size of the passive 
components, and optimization based on source and load 
specifications is required. 
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EXECUTIVE SUMMARY 


The Navy of tomorrow will require a robust and 
reconfigurable power system capable of supplying power not 
only to large high power propulsion loads but to growing 
combat system loads like high power radar and pulse loads 
such as rail guns and f ree-electron lasers. A critical 
component in such a system is the phase controlled 
rectifier. As such, the issues associated with the 
inclusion of a power electronics rectifier need to be 
addressed. These issues include input Alternating Current 
(AC) interface requirements, the output Direct Current (DC) 
load profile, and overall stability in the presence of non¬ 
linear loads. Understanding these issues and determining 
the means of assuring compatibility with a Navy all¬ 
electric ship is the focus of this thesis. 

The objective of this thesis was to investigate 
differences between multi-pulse phase controlled rectifiers 
in continuous conduction mode (CCM) with a two-pole 
inductor-capacitor ( LC ) output filter. A 6-pulse phase 
controlled rectifier with a two-pole LC output filter is 
shown in Figure 1. This circuit is the basic building block 
module for higher pulse count systems as it is easily 
stackable. 
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Figure 1. Circuit diagram of a controlled 

three-phase rectifier with a two-pole 
LC output filter. From [1]. 


The key differences that were evaluated between these 
multi-pulse controlled rectifiers were the critical 
inductance L crit (the inductance which is required to maintain 
the rectifier in CCM), the peak-to-peak ripple voltage as a 
percentage of the average output voltage 0// °^c( P p/av g ) (also 

commonly known as the DC output voltage ripple), and the 
total harmonic distortion in phase A of the line current 

%THD i . 

The determination of L crit was made by applying standard 
circuit analysis techniques, first to the 6-pulse 
controlled rectifier circuit of Figure 1, and then by 
extending the analysis to a p-pulse controlled rectifier. A 
plot of L crit versus the firing angle a is shown for a 6-, 
12-, 18-, 24-, and 96-pulse controlled rectifier in Figure 
2. From Figure 2, it is clear that L crit is lower for higher 
pulse count controlled rectifiers than for lower pulse 


xx 





is reduced 


count controlled rectifiers; the reduction in L crit 
by successively increasing pulse count. 



Figure 2. Critical inductance plot for controlled 

p-pulse rectifiers with unity resistance 
operating at 60 Hz. 

The determination of 0 /oV C ( pp/avg ) was a l so made by applying 

standard circuit analysis techniques to the 6-pulse 
controlled rectifier circuit of Figure 1 and extrapolating 
the results to a p-pulse controlled rectifier. However, 
slight approximations were introduced in the impedance of 
the load and the determination of the capacitor peak-to- 
peak ripple voltage. A plot of 0// °^c( P p/m’g) versus the firing 
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angle a is shown for a 6-, 12-, 18-, 24-, and 96-pulse 
controlled rectifier in Figure 3. From Figure 3, it is 
clear that °^c( P p/avg) is lower for higher pulse count 
controlled rectifiers than for lower pulse count controlled 
rectifiers; the reduction in °^c(ppiav g ) is dramatically 
reduced by successively increasing pulse count. 



Figure 3. Capacitor peak-to-peak ripple voltage as a 

percentage of the average output voltage for 
a 6-, 12-, 18-, 24-, and 96-pulse controlled 
rectifier with a two-pole LC output filter. 


The calculation of %THD i was performed via simulation 
using Simulink® on four p-pulse controlled rectifiers with 
the same predefined passive components. A plot of %THD i 
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versus the firing angle a for the simulated 6-, 12-, 18-, 

and 24-pulse controlled rectifiers is shown in Figure 4. 
From Figure 4, it is clear that %THD j is lower for higher 
pulse count controlled rectifiers than for lower pulse 
count controlled rectifiers; the reduction in %THD { is 
greatly reduced by successively increasing pulse count. 



versus a for p-pulse controlled rectifiers 
with a two-pole LC output filter. 
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All of these results led to the main conclusion that 


increasing the pulse count is beneficial in that it reduces 
L crit , 0 /°V C \ppiavg)' anc ^ %THD i . However, in reality, tradeoffs 

must be considered in terms of pulse count versus 
complexity and size of the power electronic and passive 
components. 


LIST OF REFERENCES 


[1] R. W. Ashton, Notes for EC3150 (Solid State Power 

Conversion), Naval Postgraduate School, unpublished, 
2012 . 


XXIV 



ACKNOWLEDGMENTS 


Without a doubt, this thesis is one of the most 
difficult undertakings I have ever accomplished in my life. 
While the work was my own, I could never have hoped to have 
completed it without the help of many people. First, I 
would like to thank my thesis advisor. Dr. Robert Ashton. 
He presented me with a thesis topic the likes of which I 
had worked on as an undergraduate student, but never fully 
understood nor could get to completely function properly. 
Undaunted, I accepted this project, inspired and hungry to 
discover and learn every detail about it that I could in my 
short time at Naval Postgraduate School (even so, I am sure 
I only scratched the surface) . His help in finding many of 
the errors I had made in the countless derivations and 
equations presented in this thesis was invaluable. Also, I 
would like to thank Mr. Andrew Parker for his role as a 
second reader of this work. 

I would like to thank my wife, Carolyn, and my 
children, Anthony and Ashley, for understanding and 
enduring the sacrifice this project has caused. 
Additionally, I always appreciated the encouragement and 
pride in me from my mother- and father-in-law. Sue and Tom 
Scaffidi (God rest his soul), and my mom, Eugenia Cray. 

I could not have done this without the friendship and 
camaraderie from some of the fine people I have met here 
who I hope to stay in touch with and remain friends with 
for the rest of my life, especially John Krbec and Karma 
and Aida Simons and Ryan, Nikki, and Gwyneth Kelly. I also 
wish to express my thanks and extend fair winds and 


xxv 



following seas to Ben and Michelle Pimentel, Tim Policar, 
Will Fleming, Jennifer Sanders, Kaylene and Brandon Carter, 
Ashley and Carson McAbbe, and Alex, Julie, and Nathan Rios. 

Finally, I would like to thank the Engineering Duty 
Officer (EDO) community for selecting me for lateral 
transfer and believing that I was a worthwhile investment 
in light of having already received one Navy-funded 
master's degree. The EDO community prides itself on being 
the technical leaders of the Navy and therefore requires 
that its members obtain technical master's degrees from NPS 
or the Massachusetts Institute of Technology. To you I say 
mission accomplished, and I will not let you down! 



I. 


INTRODUCTION 


A. BACKGROUND 

A multi-pulse phase controlled rectifier is a type of 
alternating current (AC) to direct current (DC) converter. 
This means that it takes an AC voltage source (typically a 
sinusoid) and changes it into a DC voltage. In reality, the 
output will be a higher frequency AC waveform with a 
greatly reduced peak-to-peak voltage that has an average DC 
value. The term multi-pulse usually applies to three-phase 
systems and has to do with the fact that there are multiple 
arch-shaped pulses corresponding to the peak regions of all 
of the different phase combinations of the source voltages 
that are seen by the rectifying circuit elements. By using 
transformers to create phase shifts in the source voltages 
and connecting these rectifiers together in different 
fashions, rectifiers with very high pulse counts can be 
made. A controlled rectifier means that the DC output 
voltage can be controlled or adjusted in value, whereas an 
uncontrolled rectifier will produce a fixed DC output. 

A recent and growing trend toward an all-electric ship 
and the development of integrated power systems has caused 
the role of these types of rectifiers to become even more 
important. Due to the ever-increasing shipboard power 
requirements to support such items as rail-guns, free- 
electron lasers, high power radars, and electromagnetic 
aircraft launch systems, the Navy has a tremendous need for 
rugged and reliable high power conversion modules in the 
megawatt range [1]. 
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In addition to the problem of how to meet all of the 
power demands caused by these new types of loads is the 
fact that the power electronics (the phase controlled 
rectifiers) supplying these loads are non-linear. Non¬ 
linear loads produce unwanted harmonic distortion as a by¬ 
product which can negatively affect all directly connected 
electrical equipment. The phase controlled rectifiers must 
be capable of keeping the produced harmonics to an 
acceptable level [2]. 

B. OBJECTIVES 

The differences between multi-pulse phase controlled 
rectifiers in continuous conduction mode (CCM) with a two- 
pole inductor-capacitor (LC) output filter are investigated 
in this research. The differences that are specifically 
investigated are how the critical inductance, DC output 
voltage ripple, and total harmonic distortion (THD) are 
affected by having differing pulse counts. 

C. RELATED/PREVIOUS WORK 

In [3] and [4], many of the equations for the critical 
inductance of a 6-pulse controlled rectifier with a two- 
pole LC output filter are developed and presented. However, 
they do not develop equations beyond those for a 6-pulse 
converter. Similarly, the discussion regarding DC output 
voltage ripple was limited to that of a 6-pulse converter 
in [3] and [4]. However, there has been much work and 
literature presented on the reduction of THD by increasing 
pulse count. In [5], a 24-pulse count rectifier topology is 
achieved from a single three-phase source using a novel 
interconnection of conventional single and three-phase 
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transformers to effectively produce four three-phase 
systems. In [6], an 18-pulse rectifier based on using a 
differential delta autotransformer to achieve the necessary 
phase shifts is presented for aerospace applications. In 
[7], a pulse doubling technique for a 36-pulse converter is 
accomplished via two paralleled 18-pulse converters using a 
delta/polygon-connected transformer. In order to double the 
number of pulses from 36 to 72, a tapped interphase 
transformer with two additional diodes is added at the 
output of the rectifier. In [8], the same authors present a 
similar pulse doubling technique in a 12-pulse converter 
which uses two paralleled 6-pulse converters connected 
through a star-wired transformer. A tapped interphase 
reactor is connected at the output to double the pulses 
from 12 to 24. An interesting design is presented in [9] by 
way of a parallel-connected controlled 18-pulse rectifier. 
Under reduced loading conditions, the firing angle 
increases, which results in higher THD and lower system 
efficiency. To combat this effect, the authors propose a 
dynamic control scheme which involves simultaneous and/or 
selective control of the firing angle and converter pulse 
count. For pulse count control, an appropriate 6-pulse 
rectifier or pair of 6-pulse rectifiers is disabled or 
enabled depending on the current loading. 

D. THESIS ORGANIZATION 

Background material on single-phase rectifiers, 
starting with uncontrolled diode rectifiers and concluding 
with controlled rectifiers, is provided in Chapter II. 
Critical inductance calculations are developed, first for a 
6-pulse controlled rectifier, then for a p-pulse controlled 
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rectifier, in Chapter III. The expressions and equations 
for the DC output voltage ripple, are developed in Chapter 
IV, first for a p-pulse controlled rectifier with a purely 
resistive load, then for a p-pulse controlled rectifier 
with a two-pole LC output filter. The general concepts of 
harmonics and THD are presented in Chapter V, as well as 
how they can be reduced by increasing pulse count. 
Simulink® simulations of a 6-, 12-, 18-, and 24-pulse 
controlled rectifier with a two-pole LC output filter are 
provided in Chapter VI . Conclusions and recommendations for 
future work are discussed in Chapter VII. 
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II. SINGLE-PHASE RECTIFIERS 


A. UNCONTROLLED SINGLE-PHASE RECTIFIERS 

1. Uncontrolled Single-Phase Half-Wave Rectifier 

Before being able to fully discuss what a multi-pulse 
phase controlled rectifier in CCM with a two-pole LC output 
filter is, it is useful to begin with a far simpler 
rectifier circuit commonly known as a half-wave rectifier. 
One such half-wave rectifier circuit is illustrated in 
Figure 1, where a sinusoidal voltage source is applied to a 
resistive load through an ideal diode. The diode behaves 
like a check-valve, only allowing positive AC current to 
pass through to the resistor. During the negative portion 
of the source voltage waveform, the diode blocks the 
current flow, resulting in no voltage drop across the load. 
Profiles of the source voltage, load voltage, and current 
are shown in Figure 2. From Figure 2, it is clear how the 
circuit of Figure 1 earned its common name, due to half of 
the source voltage input waveform being applied to the 
resistive load [3], [4]. 
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2 . 


Uncontrolled Single-Phase Full-Wave Rectifier 


A full-wave rectifier is the next logical progression 
beyond a half-wave rectifier. A simple full-wave rectifier 
(sometimes called a bridge rectifier) with a resistive load 
is shown in Figure 3. 


<L 



Figure 3. Uncontrolled full-wave rectifier circuit 

with ideal diodes and a resistive load. From [3]. 

In this circuit, the source voltage and load resistor 
are identical to those of the half-wave rectifier case. In 
the full-wave rectifier circuit, in addition to the current 
that flows during the positive portion of the source 

voltage waveform, current will also flow during the 

negative portion due to the second pair of diodes. During 
the positive portion of the source voltage waveform, diodes 
Z), and D 4 are "on" while diodes D 1 and D i are "off ," and 
current will flow in the direction indicated by the load 

current i L . During the negative portion of the source 

voltage waveform, diodes D { and D 4 are "off" while diodes 
D 1 and D 3 are "on." Notice that the direction of i L is 
unchanged. This unidirectional current results in a 
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positive load voltage drop during both the positive and 
negative portions of the source voltage waveform. Profiles 
of the source voltage, source current, load voltage, and 
load current as described are illustrated in Figure 4 [3], 

[4] . 



Figure 4. Source voltage, source current, load 

voltage, and load current waveform profiles for an 
uncontrolled full-wave rectifier circuit. After [3]. 


B. CONTROLLED SINGLE-PHASE RECTIFIERS 

1. Controlled Single-Phase Half-Wave Rectifier 

If the diode from the circuit of Figure 1 is replaced 
with a thyristor or silicon controlled rectifier (SCR), the 
result is a controlled half-wave rectifier of Figure 5. The 
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SCR behaves similarly to a diode in that it is a one-way 
device and will block current flow in the negative 
direction. However, it will not conduct in the forward 
direction until an appropriate trigger signal has been 
applied to its gate [10]. In reality, there will be a short 
delay before the SCR turns "on" even after it has been 
adequately triggered; but for the purpose of explanation, 
the SCR in Figure 5 can be considered ideal and will turn 
"on" immediately when triggered. The angle at which the SCR 
is triggered is commonly called the firing angle a . If the 
SCR in Figure 5 is triggered at time t l , which is equal to 

a I co (where co is the radian frequency), the resistor 
voltage and current will be as shown in Figure 6 [3], [4] . 



Figure 5. Controlled half-wave rectifier circuit 

with an ideal SCR and a resistive load. From [3]. 
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1 1 T/2 T t 

Figure 6 . Source voltage, load voltage, and load 

current waveform profiles for a controlled half-wave 
rectifier circuit triggered at time t l . From [3] . 

2. Controlled Single-Phase Full-Wave Rectifier 

By replacing the diodes in Figure 3 with SCRs, a 
controlled full-wave rectifier is created as shown in 
Figure 7. SCRs are gated in pairs with the firing angle a 
again corresponding to time t l =ala>. During the positive 
portion of the source voltage waveform, all SCRs are 
initially "off" until SCRi and SCR 4 are gated "on" while 
SCR 2 and SCR 3 remain "off" due to being reverse biased. 
During the negative portion of the source voltage waveform, 
the opposite situation occurs and all SCRs are "off" until 
SCR 2 and SCR 3 are gated "on" while SCRi and SCR 4 remain 
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reverse biased and "off." The resulting current yields zero 
voltage drop across the load resistor during the time the 
SCRs are "off" and a positive load voltage drop during both 
the positive and negative portions of the source voltage 
waveform when the appropriate SCRs have been gated "on." 
Profiles of the source voltage, source current, load 
voltage, and load current as just described are illustrated 
in Figure 8 [3] , [4 ] . 


*■ 



Figure 7. Controlled full-wave rectifier circuit 

with ideal SCRs and a resistive load. From [3]. 
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Figure 8. Source voltage, source current, load 

voltage, and load current waveform profiles for a 
controlled full-wave rectifier circuit triggered at 
time t x . t 2 = t l +T/2. After [3]. 
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III. CRITICAL INDUCTANCE 


A controlled three-phase rectifier with a two-pole LC 
output filter is depicted in the bridge configuration of 
Figure 9. In this configuration, the number of pulses is 
twice the number of phases, hence, three-phase becomes 
synonymous with 6-pulse [11] . One very useful circuit 
element parameter for these types of rectifiers is the 
value of the inductor necessary to maintain the circuit 
operating in COM. This particular value of inductance is 
called the critical inductance L crit . For nominal circuit 
operation, the instantaneous inductor current is required 
to always be positive, which is the definition of COM, and 
the capacitor is required to be very large such that the 
instantaneous capacitor voltage is approximately equal to 
the average capacitor voltage v c (t) = V c( } [3], [4]. While it 

does require lengthy equations and calculations, we can 
solve for L crit , as will be shown in the continuing sections. 





Figure 9. Circuit diagram of a controlled 

three-phase rectifier with a two-pole LC output 

filter. From [4]. 
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A. 


AVERAGE LOAD VOLTAGE CALCULATION 


To aid in understanding the relationship between the 
voltage sources for a three-phase circuit such as the one 
shown in Figure 9, the three line-to-line source voltages 
are plotted in Figure 10. While each voltage waveform 
possesses the same magnitude, they are offset from each 
other by 120°. 



Figure 10. Line-to-line source voltage waveforms 

V ab , V bc , and V ca with unity magnitudes for a 
hypothetical three-phase circuit. 

In addition to the three line-to-line source voltages 

as indicated by the polarities in Figure 9, there are three 

more line-to-line voltages that exist by looking at them 
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from the opposite polarities. These additional voltage 
waveforms, along with the originals, are plotted in Figure 
11 . 



Figure 11. Line-to-line voltage waveforms 

V cb' V ab ' V ac' V bC V ba ' and V ca wlth Unity 

magnitudes for a hypothetical three-phase circuit. 

From Figure 9 and [3] and [4], the average capacitor 
voltage V c(avg) is given by 

/- a+l20° 

V c ( a Vg )=y- J V m smydy , (1) 
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where V m is the peak value of the line-to-line voltage. The 
six in the numerator is due to the fact that there are six 
voltage pulses over one period for a three-phase rectifier, 
as illustrated in Figure 11. The integration limits are 
solved for by first dividing one period of 360° by six 
pulses to obtain a pulse-width of 60°. Dividing this pulse- 
width by two and adding and subtracting it with 90° yields 
the upper and lower limit firing angle offsets for the 
upper and lower integration limits. By integrating (1), we 
get 

V 

V C(avg ) 

Evaluating (2), we obtain 

3V r 

V c , ) =-—[cos(or +120°) - cos(a + 60°)] , (3) 

n 

which, by applying the trigonometric identity, 

cos(a + b) - cos(c/)cos(^)-sin(a)sin(&) , becomes 


31/ 


n 


-cos/ 


V=f7+t S0° 


(2: 


yj [ cos(a) cos(120°) - sin(a)sin(120°) 

m J 

n |-[cos(a)cos(60°)-sin(a)sin(60°)] 


which ultimately simplifies to 


V, 


C(avg) 


3VL 


K 


-cos a 


(5) 


B. PROCEDURE TO SOLVE FOR THE BOUNDARY ANGLE a b AS THE 
TRANSITION POINT BETWEEN SUB-MODES OF CCM 

In the circuit of Figure 9 there exists a firing angle 
a h such that the instantaneous source voltage is equal to 

the average capacitor voltage |v s (a)| =V C( ^ . Therefore, for an 
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as v. 


instantaneous source voltage v s (7) given as v s (f) = V m sin(<yt) and 

an average capacitor voltage previously solved for in (5) , 
we have the equality 

3V 

V m sin(fei/,) = —- cos a h ( 6) 

n 

By substituting <^=< 2 ^ + 60°, we have 


3VL 


V m sin(« fo + 60°) = —— cosa ( 


n 


(7; 


By applying the trigonometric identity 

sin(a+b) = sin(a)cos(Z?) + cos(a)sin(Z?), we produce 


3V„ 


V m [sin ( a b ) cos (60°) + cos ( a b ) sin (60°) ] = —— cos a b , 


which can be rewritten as 


— sin a, + V3 cos a, 
2 L 


= —cosa b 
n 


Grouping similar terms yields 


sin a h = 


f—- V3 


cos a, , 


: i o: 


which can then be reduced, leaving 


tan a h = 


f—- V3 


( 11 ) 


which can be solved for a h to two decimal places as [3], [4] 


a, = 10.08° . 


12 ' 
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C. PROCEDURE TO SOLVE FOR THE CRITICAL INDUCTANCE FOR A 

FIRING ANGLE > 10.08° 

With a continuous inductor current and for firing 

angles greater than 10.08°, the absolute value of the 
instantaneous source voltage is greater than the average 
capacitor voltage at the point of commutation, given by 
a + 6ti° = cot l . Therefore, the slope of the inductor current at 
the point of commutation is always greater than zero, and 
it can be said that the voltage source is always 

electrically connected to the load via the SCR bridge 

network. In Figure 12, the waveform profiles of the 
instantaneous load voltage, average load voltage, and 
instantaneous inductor current for the circuit of Figure 9 
are provided. In order to establish a well-defined boundary 
condition for the proceeding calculations, we let L = L crjl so 
that i L (t l ) = i L (t 2 ) = 0 every one-sixth period of the source 

voltage, resulting in operation on the boundary between 
continuous and discontinuous modes (also known as being 
barely continuous), as shown in Figure 12. 
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Figure 12. Load voltage, average load voltage 


and inductor current waveform profiles for a 
controlled three-phase rectifier with a two-pole LC 
output filter with a firing angle greater than 10.08°. 

After [3]. 

From Kirchhoff's Voltage Law (KVL), we have 



(13) 


where v 5 (V) is the source voltage given by v s (?) = V m sin(a>f) . 


Substituting (5) for V C{avg) and L for v x (t) , we get 



di, (t) 3V 

V m sm(®/) = L -— + —— cos a 

dt n 


(14) 


Rearranging (14), we get 



dt L\ n 


y • 


(15) 
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To solve for i L (t ), it is necessary to set up the integral of 
(15), leading to 


V r 


sin(<z>r)-cos a 

n 


dt , 


:i6; 


which, when integrated, results in 


V 

i L (t) = -r 


-cos (cot) 3tcosa 


co 


n 


Jf, * 


(17) 


Upon evaluating (17) and factoring out 1/co, we get 


V 

i L (t) = 

coL 


-cos (cot)- 


3cotcosa 


n 


3 cot, cos a 


— COS( GJt\) — —7 

v n J 


(18) 


Rearranging (18), we get 


V 

hi?)-= 

coL 


f 1 \ 


cos(<zw\) - cos (cot) + j(oot\ - cot ) cos a 


:i9; 


By substituting cot=a + 6 0° in (19), we get 


V 

coL 


cos(a + 60°)-cos(<ztf) + 


i 3 V an n ^ 
+- cot 


\n j 


180 3 


cos a 


(20: 


It is now necessary to calculate the average load current 
hvg ' given by 


avg rj. 


1 T 

-\h(t)dt 


(2i; 


Substituting (20) into (21), we get 


/ - 6V »- 

flV * 2 ncoL 


a+ 120 ° 


a+6 0 ° 


cos(a + 60°) - cos (cot) + — 

n 


3 f an n ^ 
■ +- cot 


180 3 


cos a 


d(cot) . (22) 
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Integrating and evaluating the integral given by (22) in 
one step would be very cumbersome, so it is broken up into 
three parts. Let 


a+120° 

X x = cos(a + 60°) J d(cot), 

a+ 60 ° 

a+ 120 ° 

X 2 = J [-cos(fflf)]rf(fflf), 

a+ 60 ° 


and 


X 3 = 


3 cos a 


n 


a+ 120 ° 


a+ 60 ° 


an n 

-+- cot 

180 3 


d(cot) 


Starting with X l and integrating (23), we get 


X x = [<wfcos(eir + 60 0 )] 


cot=a+ 120 ° 
cot=a+6 0 ° 


which, when evaluated, becomes 


X x =ycos(a + 60°) . 

By applying the trigonometric identity 

cos(a + b) - cos(a)cos(^)-sin(a)sin(&) to (27), we get 


X x =y[cos(a)cos(60°)-sin(a)sin(60°)] , 


which can be rewritten as 


X, = 


K 


cos(a)-V3sin(a) 


Continuing with X, and integrating (24), we get 


r • / ^-\cot=a+ 120° 

X 2 = [-sm(<»f)] , 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


( 30 ) 


which, when evaluated, becomes 
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X 2 =sin(a + 60°)-sin(a + 120°) . 


'31' 


By applying the trigonometric identity 

sin(a + b) = sin(a)cos(£) + cos(a)sin(&) to (31), we get 

X 2 = sin(a) cos(60°) + cos(a)sin(60°) 
-[sin(a)cos(120°) + cos(< 2 )sin( 120 °)] 


'32' 


which simplifies to 


X 2 = sin a 


' 33) 


Finishing with X , and integrating (25), we get 


*3 = 


3 cos or 

V 7t J 


f an n\ [cot) 


180 + 3 


cot-- 


-\cot=a+\20 o 


-io)t=a+60° 


'34' 


which, when evaluated, becomes 


*3 = 


3 cos a 
y 7i j 


f OK 7t\( 7t\ 1 | ( OK 2 K 

+ — 


180 3 


n 

v 3 j 


v 


180 + ’ 3 


OK K 
■ + — 


180 3 


:35; 


Simplifying (35), we get the result 


*3=~ 


^■cosa 


:36; 


Inserting the resultant expressions from (29), (33), and 

(36) back into (22), we have 


6V m k 
m ' g 2 kcoL \ 6 


cos(a)-V3sin(a) 


+ sina-- 


^■cosa 


'37' 


which, when simplified, produces 

3V„, sin a 


I avg 


kcoL 


\_kS' 


os: 
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Because 

we have already 

solved 

f ° r V C{mg) in 

(5) , we have 

another 

expression for I avg 

using 

Ohm's Law, as 



Kvg =—^-cosa. (3 9) 

7rK 


By substituting (39) into (38), we can solve for L crjt as [3], 
[4] 


^crit 


Rtma 


co 


n 


■S 


(40) 


D. PROCEDURE TO SOLVE FOR THE CRITICAL INDUCTANCE FOR A 

FIRING ANGLE < 10.08° 

With a continuous inductor current and for firing 
angles less than 10.08°, the absolute value of the 
instantaneous source voltage is less than the average 
capacitor voltage at the point of commutation. Therefore, 
the slope of the inductor current at the point of 
commutation is always less than zero. Unlike the previous 
sub-mode in which the minimum inductor current occurred at 
the firing point, for this sub-mode, the minimum inductor 
current occurs at an angle 6 3 = cot i (where 6 = # 3 - 60 ° referenced 
to a ), as shown in Figure 13. 
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Figure 13. Load voltage, average load voltage, and 
inductor current waveform profiles for a 
controlled three-phase rectifier with a two-pole LC 
output filter with a firing angle less than 10.08°. 

After [3]. 


Once again, we can say that the voltage source is 
always electrically connected to the load via the SCR 
bridge network, and the conducting SCRs are able to ride- 
through the minimum current point at t 3 . At t 3 , the average 
capacitor voltage is equal to the instantaneous source 
voltage and the slope of the inductor current is zero. 
Therefore, for an instantaneous source voltage v s (?) given 

as v,(f) = V m sin(fflf 3 ) and an average capacitor voltage previously 
solved for in (5), we have the equality 


V m sin(<af 3 ) 


3VL 


n 


-cos a 


(41) 


, we have 


24 


By substituting ®/ 3 =6 , + 60° 









V„sm(0+6O°) 


(42) 


3VL 


K 


-cosa , 


which can be solved for 9 as 


9 = sin 1 


— cos a 
\n j 


-60 c 


(43) 


For t>t 3 , the instantaneous source voltage is greater than 
the average capacitor voltage, and the inductor current 
begins to increase. By triggering the SCRs at a b as 
previously determined in (12), an additional verification 
can be made at the boundary between sub-modes. For the 

instantaneous source voltage v s ( 7 ) = V m sin(taf 3 ) and an average 
capacitor voltage as previously solved for in (5) (but with 
a b now substituted for a ) , we have the equality 


By substituti 


which can be 


V,n sin (cot ,) = 
ng cot 3 = 0 h + 60° , we 
14,801(3,+60°) 
solved for 0 b as 


3VL 


cos a, . 


n 

have 



cos a b , 


(44) 


(45) 


0 h - sin 1 


^ 3 

— COS (2, 

\n 




-60° . 


(46) 


At a h = 10.08°, 6 h = 10.08°. In order to once again establish a 

well-defined boundary condition for the subsequent 
calculations, we let L = L crit so that i L (t 3 ) = 0 every half-period 
of each phase of the source voltage, again resulting in 
operation on the boundary between continuous and 
discontinuous modes, as illustrated in Figure 13. We begin 
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once again with the KVL equation as given in (13), with the 
same expressions for v s (?) , v v (?), and V C{avg) r working through 

identical equations given by (14) and (15), and arrive at 


hit) 


V 




3 

-cos a 

n 



(47) 


After integration, we have 


and after 
result 


hit) 


evaluating 


V 

_ m 

" L 
(48) 


-cos (cot) 3t cos a 

I 

co n 

-ir 3 

and factoring out 


1 I CO , 


(48) 

we get the 


y 

h(t) = ^T 
coL 


-cos (cot)- 


3cotcosa 


K 




3cot 2 cosa 
n J 


■49: 


Rearranging (49), we get 


y 

4(0 —y- 

coL 


f 3 i 

cos(< r 6 » L 3 )-cos(ft>/ 1 )+ — (<j# 3 -cot)cosa 

Vx ) 


:5o: 


By substituting cot^=0 + 6O° into (50), we get 


V 

hit) = ~ 

coL 


cos (0 + 60°) - cos(a »0 + 


7-2 1 




On n 


- H- cot 


180 3 


cos a 


:5i: 


Once again, we calculate the average load current using 
(21). By substituting (51) into (21), we get 


I 

avg 


6V 

_ m 

2 nco L 


a+ 120 ° 


a+6 0 ° 


cos ( 6 * + 60°) - cos (cot) + — 

n 


C On n 

-+- cot 

v 180 3 


cosa 


d(cot ) 


(52) 


Integrating and 

evaluating 

the 

integral given by ( 

[ 52) in 

one step 

would 

once again 

be 

very cumbersome, so 

it is 

separated 

into three parts. 

Let 
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a+ 120° 

X x = cos(<9 + 60°) J d(cot ), 

a+ 60° 


X 1 =X 1 from (24), and 


*3 = 


3 cos a 


n 


a+ 120 ° 


a+ 60° 


On n 
- +- cot 


180 3 


d(cot ) 


( 53 ) 


(54) 


Starting with X l and integrating (53) , we get 


[ / ^ ~\cot=a+\2W 

cot cos(0 + 6O')1 

V /Jfflr=a+60° 


which, when evaluated, becomes 


X t =yCos(0 + 6O°) 


(55) 


(56) 


Since X 2 is identical to the previous X , , the result is 
provided from (33) . Finishing with X 3 and integrating (54), 
we get 

ox=a+ 120° 

, (57) 

art=a+ 60° 


*3 = 


^3cos a^ 

\ n y 


On n 


180 3 


cot- 


(ax)' 


which, when evaluated, becomes 


*3 = 


3 cos a 


V n ) 


On n 
180 3~y 




j 


1 I ( an In 

vl80 T 


y r \ 2 

' ' an n ' 

v i8o y 


1 j 


Simplifying (58), 


we get 


*3 = 


^ n "3 


Vl 80 y 


(#-a-30°)cosa . 


(59) 


Inserting the resultant expressions from (56), (33), and 

(59) back into (52), we have 
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( 60 ) 


n ( \ 

/ =—— —cos (0 + 60°) +sin or + (0-ar -30°) cos a 

avg ttcoL [_ 3 V 7 U 80 j V 7 

Using (39) once again, we can substitute it into (60) and 
solve for L crit as 

L =——— — cos(# + 60°) + sina + —— (6-a -30°)cosa (61) 

cocosa [ 3 ^180 y J . 

Although 9 is a function of a as given by (43) , 
substituting it into (61) would only serve to make it more 
unwieldy, so we leave it in its present form [3], [4]. 

E. SUMMARY OF CALCULATIONS FOR L :nt 

To summarize, we have the following equations for L crjt : 



R tan a n V3 

® l 6 , 


for a >10.08° 


(62) 



From Figure 14, the critical inductance for the 6-pulse 
rectifier circuit of Figure 9 can be seen, assuming a unity 
resistance and a frequency of 60 Hz. An important 
conclusion is that as a increases, L crit increases at an 
increasing rate. This is the key reason that circuit 
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designers of controlled rectifiers desire to keep a 
relatively small—the size, weight, and cost of the inductor 
cannot be ignored. 



0 10 20 30 40 50 60 70 80 90 

a (degrees) 

Figure 14. Critical inductance plot for a controlled 
6-pulse rectifier with unity resistance 
operating at 60 Hz. 

F. CRITICAL INDUCTANCE FOR A CONTROLLED P-PULSE RECTIFIER 

One way to increase the pulse count of this type of 
rectifier is by stacking 6-pulse rectifiers in series—for 
instance, two 6-pulse rectifiers could be connected 

together to make one 12-pulse rectifier. While this has 
several benefits, an important one is that it decreases the 
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value of L crit . To prove this, one could certainly go back and 
tediously derive the necessary equations, solve for L crit , and 
plot it. Indeed, that is exactly what was done for the 12-, 
18-, and 24-pulse rectifier circuit configurations, the 
workings of which have been spared from this thesis. The 
purpose of this labor was to uncover a possible pattern. 


This 

work 

was not done 

in vain; 

a clear 

pattern exists. 

There 

are 

equations to 

calculate 

L crit f ° r 

a controlled p- 


pulse rectifier. These equations are derived in the same 
manner as before, only this time with some clever 
substitutions. First, let p equal the pulse count. Then, 

let the lower integration limit 0 L be expressed in degrees 
as 


0 L = 90 ° - 


180 ° 

P J 


65! 


and the upper integration limit d u , also in degrees, be 

180 ° 


e v =9o°+ 


V P J 


66 ! 


thereby, making both functions of pulse count. Next, let 
A = cos 0 l , B = sin 9 l , C = cos0 v , and D = sind u ; therefore, C = -A 

and D-B . 


1. Average Load Voltage Calculation 

With our new substitutions, the expression for the 
average load voltage from (1) now becomes 

a+6jj 

v c(a Vg )=^~ j K.sin rdr' (67) 

lK A 
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which, when integrated, yields 


V c( avg ) =-^ L [co^a + 0 u )-co^a + e L )] . 
in 


68 ' 


By applying cos(a + b) - cos(tf)cos(7>)-sin(a)sin(7>) to (68), we get 




C(avg) 


- ^ m [Ccos(a) -Dsin(a) - Acos(a) + fisin(a)], 


2 n 


691 


which simplifies to 


v aa V g) = —^—^cosa 
n 


;7o; 


2. Procedure to Solve for the Boundary Angle a b as 
the Transition Point Between Sub-modes of CCM 

By beginning with (6), but this time substituting 
tot l =a h + 0 L and our new expression for average capacitor 
voltage from (70), we have 

V„,sin( a h + 8,)= cos a h . (71) 

n 

By applying the trigonometric identity 

sin(a + b) = sin(c/)cos(&)+ cos(a)sin(&) , we produce 

y m [sin(a i )cos( 6 ' i ) + cos(< 2 / ,)sin( 6 ' i )] = ^^ m -cosa b , (72) 

which can be rewritten as 


Asin(a h ) + B cos(a h ) 


Ap 

— cos a h . 
n 


Grouping similar terms, yields 


Asin(a Z; ) 


fhi-A 

V n j 


cos a b , 


(73) 


(74) 
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which can then be reduced, leaving 


, \ (p B 

tan (a b ) = -- 

V 7t A 


Therefore, we have 


_i [ p B 

a, = tan- 

l n A 


3. Procedure to Solve for the Critical Inductance 
for a Firing Angle > a h 

Starting with (13), we use the same original 
substitutions, with the exception of our new expression for 
V C{avg) , and obtain 

dl yl n 1/ 

V m sin(cot) = L —- +- -cos a, (77) 

dt n 

which, after we solve for di L /dt, integrate, and simplify as 
before, we get 


y r f Ap'] 

h (0 = — 1 — cos (o)t x ) - cos (cot) + — ( cot { - cot ) cos a 
coL \ n ) 


By substituting cot l =a + 8 L in (78), we get 


i L (t) = —— cos(a + 0 L ) - cos(cot) + 


Ap\( an 0 L n 


n A180 180 


+ —- cot cosa 


Calculating I by substituting (79) into (21), we get 


cos(a + 0 L )-cos(fflf) 

hv g = ^ A | (Ap V an 0 L n \ d(cot) 

2ncoL * + — -+ -^- cot cos a 

L l n 180 180 
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Breaking up the integral into similar parts as before, we 
have: 


Xj=cos (a + 0 L ) | d(cot), 


“ ' W U 

X 2 = | [-cos(<z>f)]<i(<ztf), 


v Apcosa f" (an G.n ] 

X,=— - -+ —- cot d(cot) 

n J ll 80 180 Y 


Starting with X l and integrating (81), we get 


X l = [<jtfcos(a + # L ) 


cot=a+9u 


which, when evaluated, becomes 


x i = — (@u -9 L ) C0S (a + & L ) . 

V ioU y 


We apply cos(a + b) = cos(a)cos(&)-sin(c<)sin(&) to (85), yielding 


*i = — (^-^)[cos(a)cos(^)-sin(a)sin(^ L )], 
V ioU 


which can be rewritten as 


X t = — [~Acos(a)-.6sin(a)l 

VP J 


Continuing with X 2 and integrating (82), we get 


*.=[-“ <«)& 


which, when evaluated, becomes 
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X 2 = sin(a + 0 L )~ sin(« + 0 V ) . 


: 89 ; 


By applying sin(c< + &) = sin(c/)cos(£) + cos(c/)sin(&) to (89), we get 

X 2 = sin(a)cos(6 , L ) + cos(< 2 )sm(# L ) 

- [sin(< 2 ) cos^) + cos(a) sin^)] 


90 : 


which simplifies to 


I, = 2Asin<2 


91' 


Finishing with X 3 and integrating (83), we get 


X 3 = 


Ap cos a 

V n ) 


an 6,n 
■+ L 


180 180 


^ M 

cot --—— 

J 2 


cot=a+0,, 


cot=a+Q, 


92' 


which, when evaluated, becomes 


X 3 = 


Ap cos a 

\ n j 


( an 0 T n ^ 

\{O V -0 L )\ 

f n ^ 

-+ — 


OO 

o 

00 

o 

1 V u L / | 

U8oJ 


1 I (an 9„n\ ( an 9,n 2 2 

■+ L 


180 + 180 


180 180 


93: 


Simplifying (93), the result is, 

_ 2 An cos a 

2^-3 — 


94' 


Inserting the resultant expressions from (87), (91), and 

(94) back into (80), we have 


I - pV "< 
avg IncoL 


/ 2n^ 

l P o 


[A cos a -5 sin a] + 2Asina - 


2 An cos a 


95 ; 


which, when simplified, produces 

V„ sin a 


avg ncoL 


( Ap-nB ) . 


96; 
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in (7 0), we have 


Because we have already solved for V 


C(«vg) 


another expression for I using Ohm's Law, as 


, ApV m 
J a Vg = — -f-cosa . 

71K 


97 ; 


By substituting (97) into (96), we can solve for L . as 


L V = 

crit 


Rian a 


f 


co 


1 - 


7rB 

Ap 


98' 


4. Procedure to Solve for the Critical Inductance 
for a Firing Angle < a b 

Starting with (41), we substitute our new expressions 
for V C(avg) and cot 3 to obtain 


V m sin (6 + 0 L ) = ApV>" cos a 

K 


99; 


which can be solved for 0 as 


6 = sin 


Ap 


-cos a 


-dr 


:ioo: 


Once again, we can make an additional verification at the 
boundary between sub-modes. For the instantaneous source 
voltage v s (t) = V m sin(<5# 3 ) and an average capacitor voltage as 
previously solved for in (70), we have the equality 


K,sin (fi# 3 ) = 


A PV m 


n 


cos a b . 


101 ' 


By substituting ^ 3 = 6 *,, + ^ into (101), we have 


V m sin (d h + d, ) = cos a b , 
n 


102' 


which can be solved for 6, as 
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To calculate the average load current, we again refer to 
(21). By Substituting (108) into (21), we get 

T PV„, t n n \ Ap ( On 6,71 \ 1 

I = — cos(6 + 6, )-cos(&)r) H -- +~ 0Jt cos a dicot) (109) 

flV? IkojL J di l v LJ n {ISO 180 J . 

Because of the complexity of (109), it is again broken into 
three parts. Let 
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( 110 ) 


(X-\- 0 r/ 


X x = cos(# + # L ) | d(cot) , 


a+6. 


= X 0 from (82), and 


Ap cos a a r' : ( Ok 0,k 
X, = —- -+ -^- cot 

k JAm iso . 


3 


d(cot) 


Starting with X x and integrating (110), we obtain 


X x = '(cot cos(0+0^ 


cot=a+G u 

cot=a+0 L 


which, when evaluated, becomes 


X, 




V 180y 


(0 u ~0 L )cO&(0 + 0 L ) f 


which can be rewritten as 


*i = 


f In'' 

P , 


cos 


(0 + 0l) 


Because X, is the same as the previous X, , its 
given in (91). Finishing with X 3 and integrating 
obtain 


*3 = 


Ap cos a 


V x J 


On 0,n 


180 180 


(ot- 


(cot) 


cot=a+0,, 


cot=a+0. 


which, when evaluated, becomes 

I ^ On 0,n^ 

f Ap cos a ^ 


X 3 = 


\ n ) 


180 + 180 




^ n ^ 


v. 180 j 


1 \( an 0,,n\ ( an 0,n ^ 

-+ L 


180 180 


180 180 


By simplifying (116), the result is 


(111) 

( 112 ) 

(113) 

(114) 

result is 
(111), we 

(115) 


(116) 
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On 

an 

y 

n 

^ 180 

180 

PJ 


By inserting the results from (114), (91), and (117) back 

into (109), we have 


, pV„ („ x . . A On an n 

/ =—— — cos [6+0, +Asina + A- cosa 

* na>L\p) v ^ 180 180 p) 


By using (97), we can substitute it into (118) and solve 


for L crit as 


, R I ft l / „ „ N | On an n) 

L.= - — cos(0 + 0, ) + Asina + A-cos a (119) 

c coAcosa [pi v ’ 1 180 180 p) 


Although 6 is a function of a as given by (100), 
substituting it into (119) would only serve to make it more 
unwieldy, so we leave it in its present form. 

5. Summary of Calculations for L crjt 

To summarize, we have the following equations for L . : 


R tan a ( nB 

L cri,= - 1-— for a > a h 

co l Ap J 


T R (n ] ,. . >. , . a On an n\ 

L.= - — cos (0 + 0, +Asina + A-cos a (121 

c " coAcosa {p J v L) U 80 180 Pj 


for a<a h , where 


if P B 


a h = tan-, 

\n A 


a ■ -1 ( A P ) a 
9 -sin —cosa -0, , 

l 71 


( 123 ) 
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(124) 


A = cos 0 L , 

and 

B = sin 0 L . (125) 

From Figure 13, critical inductances for a 6-, 12-, 18-, 

24-, and 96-pulse rectifier can be seen, assuming a unity 
resistance and a frequency of 60 Hz. As before, as a 
increases for each rectifier, L crit increases at an increasing 
rate. Another important conclusion from Figure 13 is that 
as pulse count increases, L crit decreases at a decreasing 
rate. This is most obvious when comparing a 6-pulse 
rectifier to a 12-pulse rectifier, where there is a large 
reduction in L crit , and then continuing to an 18-pulse 
rectifier, where the percent difference in L crit is less 
substantial. To clearly illustrate the concept of 
diminishing returns, the jump from a 24-pulse rectifier to 
a 96-pulse rectifier was purposely plotted to show just how 
little reduction in L crit is obtained by dramatically 
increasing pulse count. Finally, as pulse count increases, 
a b decreases, making the equation for L crit for a <a b less 
relevant, especially for high pulse count rectifiers. 
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x 10" 



Figure 15. Critical inductance plot for controlled 
p-pulse rectifiers with unity resistance 
operating at 60 Hz. 


Critical inductance values for the L crit plots of Figure 
15 are provided in Table 1, starting at a firing angle of 
0° and ending at 80°, in 10° increments. A firing angle of 
90° would result in an infinite critical inductance, so the 
upper bound is not inclusive of 90°. 
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Table 1. Critical inductance (Henry) for selected firing 
angles for the controlled p-pulse rectifiers of 

Figure 15. 



firing angle alpha 

degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

2.51E-05 

4.35E-05 

8.99E-05 

1.43E-04 

2.07E-04 

2.94E-04 

4.28E-04 

6.79E-04 

1.40E-03 

12 

3.07E-06 

1.07E-05 

2.22E-05 

3.51E-05 

5.11E-05 

7.26E-05 

1.05E-04 

1.67E-04 

3.45E-04 

18 

9.07E-07 

4.76E-06 

9.82E-06 

1.56E-05 

2.26E-05 

3.22E-05 

4.67 E-05 

7.42E-05 

1.53E-04 

24 

3.82E-07 

2.67E-06 

5.52E-06 

8.76E-06 

1.27E-05 

1.81E-05 

2.63E-05 

4.17E-05 

8.60E-05 

96 

5.96E-09 

1.67E-07 

3.45E-07 

5.47E-07 

7.95E-07 

1.13E-06 

1.64E-06 

2.60E-06 

5.37E-06 


A much more useful table would be one that shows a 
relative comparison between the absolute L crit values of 
different controlled p-pulse rectifiers, such as Table 2. 
In Table 2, the 6-pulse controlled rectifier is treated as 
the basis to which the other p-pulse controlled rectifiers 
are compared. The L crit ratios become fixed percentages beyond 
the boundary angle for each controlled rectifier. 


Table 2. Critical inductance as a percentage of 6-pulse 
critical inductance for selected firing angles for 
the controlled p-pulse rectifiers of Figure 15. 



firing angle alpha 

degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

12% 

25% 

25% 

25% 

25% 

25% 

25% 

25% 

25% 

18 

4% 

11% 

11% 

11% 

11% 

11% 

11% 

11% 

11% 

24 

2% 

6% 

6% 

6% 

6% 

6% 

6% 

6% 

6% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 
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IV. DC OUTPUT VOLTAGE RIPPLE 


The discussion of critical inductance began with the 
assumption that the instantaneous capacitor voltage was 
approximately equal to the average capacitor voltage (which 
was equal to the average load voltage). While it was useful 
then to conceptually think of the capacitor as being 
infinitely big and causing v c (t ) to be constant, we know this 
is not the case. The load voltage will have some amount of 
ripple, but it will be mitigated by the capacitor during 
its charge and discharge periods. 

A. OUTPUT VOLTAGE RIPPLE FOR A CONTROLLED RECTIFIER WITH 

A RESISTIVE LOAD 

While the concentration of this thesis is not 
specifically about phase controlled rectifiers with purely 
resistive loads, it is a good starting place to understand 
about output voltage ripple. The circuit of Figure 16 is an 
example of a controlled three-phase (6-pulse) rectifier 
with only a resistive load. However, let us consider it to 
be a controlled p-pulse rectifier for the proceeding output 
ripple voltage analysis. 
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Figure 16. Circuit diagram of a controlled three-phase 
rectifier with a resistive load. From [4]. 


From [3] and [4], for a<0 L , the instantaneous load voltage 
is always positive or zero, and its average load voltage is 
synonymous to (70) and is given by 

V = ApV,n cos a (12 6) 

n 

For a>0 L , the instantaneous load voltage becomes negative 
in parts of the cycle and is expressed as 

180 ° 

V avg =^~ J v „,smr d r' ( 127 ) 

' 2n aie L 

where the upper integration limit is replaced with 180° to 
limit the instantaneous load voltage to either zero or 
positive values. Integrating and evaluating (127), we get 

Kv g = cf^[l ■+ cos(a + d,)\ (12 8) 

2n 
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From [4], the root mean square (RMS) value of 
voltage is obtained by applying the definition 
given by 


V, 


RMS 





0 


Starting with the first case, we have 


V, 


RMS 



a+Ofj 

J sin 2 ydy 

a+0 L 


By applying the trigonometric identity skT0 = [l 
(130) and integrating, we get 


-cos 


VrmS 


IpK, 

8 7V 


1 y=a+e v 


■[2r-sin(2y)]^ 


which, when evaluated, becomes 


Vr ms = ~ ~ [ sin ( 2a + 20 u ) - sin ( 2 « + 20 l )] j 


By applying sin(a + b) = sin(c/)cos(^) + cos(a)sin(&), we have 


Vrms 


pvj 

87V 


-h sin ( 2<2 ) cos(2$ L ) + cos(2<3t) sin (20^ 

- sin ( 2a ) cos ( 26 v ) - cos ( 2a ) sin ( 20 u ) 


Next, let 
therefore. 


E = cos( 29 l ), F = sin(26* i ), G = cos(2 0 V ), and t 

H = -F and G=E. Then, (133) reduces to 


V, 


RMS 


= P ] , Fpcos(2a) 


>/2 


2 TV 


Continuing with the second case, we have 


the load 
of RMS, 

(129) 

(130) 
(20)]/2 to 

(131) 

(132) 


(133) 

= sin(26»,) ; 

(134) 
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v 

v RMS 


pK 


2 180 ° 


2 n 


J sin 2 ydy 


a+6, 


We again apply sin 2 0 — [l-cos(2#)] / 2 to this result 
integrate to obtain 


V 

v RMS 


IpK 

8 n 


_ [2/-sin(2f)] 


^= 180 ° 

y=a+0 L 


which, when evaluated, becomes 


V 

v RMS 


pK, 


8u 


ecu 2 u . / „ \ 

u - + — + sin ( 2a + 20, 

90 p v L> 


By again applying sin(a + &) = sin(ci)cos(&) + cos(a)sin(&) , we have 


Vrms 


pvj 

8 u 


„ an 2 n 

2 n - n-\ - 

90 p 


+ sin (2«) cos [20 L ) + cos (2a) sin [20 L ) 


Next, let E = cos( 20 l ) and F = sin(2# L ). By substituting 
into (138) and rearranging, we get 


_V m 2 + p ap Epsin(2a) Fpcos(2a) 
Vrms ~ V 4 360 + 4~n + 4 ~n 


Now we can compute the RMS ripple voltage V RMS _ rippIe for 
case as 


v = lv 2 -V 2 

v RMS—ripple Sj v RMS v avg 


For the first case ( a<0 L ), we substitute (134) and 
into (140) to obtain 


(135) 
and 

(136) 

(137) 

(138) 

these 

(139) 

each 

(140) 

(126) 
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V, 


RMS-ripple 


= V„I- 


1 + 


Fpcos{2a) 


2 n 


A 2 p 2 2 

- 2 —cos a , 

n 


(i4i; 


while for the second case (a>0 L ) , we substitute (139) and 
(128) into (140) to obtain 


v =V 

v RMS-ripple y m 


if 2 + p ap £psin(2a) Fpcos(2a) 




360 


4 n 


4 n 


J 


'142' 


An' 


[l + cos(« + 6 , L )]" 


Computing the peak-to-peak ripple voltage V pp - ripp i e and the 
peak-to-peak ripple voltage as a percentage of the average 
output voltage %V pp/avg are slightly more involved because 

there are actually three separate regions we must evaluate. 
The first two regions are actually sub-regions of the 
aforementioned first case, while the final region 
correlates with the second case from earlier. For the first 

region, given by a<{B u — B I ^t2 l the maximum voltage is V m , 

and the minimum voltage is V m sin (<2 + 6^) . Therefore, V pp - ripp i ei is 
their difference. 


Vpp-rippie, = K, [l - sin (a + 0 ,,)] , 


(143; 


which can be written in terms of p as 


V . , =V 

pp-ripple x m 



1 

O 

O 

00 

1-sin 

a +90 + 


l P )_ 


(144; 


The peak-to-peak ripple voltage as a percentage of the 
average output voltage is given by 


%V. 


pp/avg 


V, 


pp-ripple 


100 % . 


(145) 
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By substituting (144) and (126) into (145), we get 


1-sin a+ 90° + - 


100 % 


The second region is given by {9 V —6^)12<a<9 h , in which the 
maximum voltage is V r m sin(a + 0 i ) , while the minimum voltage is 
V m sin(ar + ^ £/ ) . Therefore, V pp _ ripp i ei for this second interval is 


pp-ripple 2 


: K [sin {a + 6,) - sin (a + d (J )] , 


which can be written in terms of p as 


f 1 QQO \ f i oAO 

y pp - ripp ie 2 =y m sin « + 90°-sin a + 90° +- 

V P J V P 


By substituting (148) and (126) into (145), we get the 
peak-to-peak ripple voltage as a percentage of the average 
output voltage 


f 1 QQO \ / 1 QQO 

sin a + 90° - sin a + 90° +- 


100 % 


(149) 


Finally, the third region is given by a>0 L . In this region, 
the maximum voltage is V' ffl sin(a + ^ L ) , and the minimum voltage 
is zero, making V pp _ ripple equal to the maximum voltage or 

y pp -n PP ,e 3 =V m sm(a + 0 L ), (15 0 ) 


which can be written in terms of p as 
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(151) 


f 


V ■ , =V 

pp-ripple 3 m 


sin 


a + 90° 


180° ^ 

P y • 


By substituting (151) and (128) into (145), we get the 
peak-to-peak ripple voltage as a percentage of the average 
output voltage 


%V. 


PP/avg 3 


f 


sin 


a+ 90°- 


180 c 

P J 


^[1 + cos (a + 0 L )] 
in 


100 % 


(152) 


V J 

The results of both cases for a 6-pulse rectifier with 
V m equal to 1.0 V and a purely resistive load are 
summarized in Figures 17 and 18. Initially, as a increases, 
the difference between the RMS load voltage and the average 
load voltage increases, making the square root of the 
difference of their squares grow, causing the RMS load 

voltage ripple to steadily increase. However, as a 
continues beyond 60°, the square root of the difference of 
their squares begins to decrease, causing the RMS load 

voltage ripple to peak and also begin to decrease. The 
peak-to-peak ripple load voltage rapidly increases until 
reaching a subtle inflection point at a- 30° and then 
increases less rapidly until reaching a maximum at a- 60°, 
where it begins to rapidly decrease. 
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0 20 40 60 80 100 120 

a (degrees) 

Figure 17. Average voltage, RMS voltage, RMS ripple 

voltage, and peak-to-peak ripple voltage for a 6-pulse 
rectifier circuit with a unity V m and a purely 

resistive load. 

The peak-to-peak ripple voltage as a percentage of the 
average output voltage increases as a increases to 60° (the 
transition at 30° is virtually undetectable) because V avg is 
fairly high and decreasing while V pp _ rjpple is low but growing 
quickly. Beyond 60°, %V pp/avg increases at a tremendous rate 
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This 


due to V decreasing at a faster rate than V pp _ ripple 

illustrates another reason why circuit designers prefer to 
keep a fairly small. 



Figure 18. Peak-to-peak ripple voltage as a 
percentage of the average output voltage for a 6-pulse 
rectifier circuit with a unity V m and a purely 

resistive load. 

In order to demonstrate how increasing pulse count 
improves the peak-to-peak ripple voltage as a percentage of 
the average output voltage for these types of rectifiers, 
%V p p/ avg has been plotted for several p-pulse rectifier 

circuits in Figure 19. The plot exhibits very similar 
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behavior to Figure 15 in that increasing pulse count yields 
better results (with rapidly diminishing returns), while 
increasing a worsens the results. While not shown in Figure 
19, further increasing a beyond 0 L (or much beyond 30°, for 
that matter) is purely academic as %V pp/avg very rapidly 

reaches exorbitant amounts. This also illustrates why a 
rectifier with a purely resistive load is an unlikely 
candidate for a rectifier which necessitates a low %V pp/avg as 

it will require a high pulse count (increased hardware 
complexity, cost, and space). 



Figure 19. Peak-to-peak ripple voltage as a percentage 
of the average output voltage for p-pulse 
rectifier circuits with unity V and 

m 

purely resistive loads. 
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To help in understanding some of the items associated 
with angles for the p-pulse rectifiers with purely 
resistive loads, a quick reference summary is provided in 
Table 3. For instance, for a 12-pulse rectifier with a 
purely resistive load, its first pulse starts at 75° and 
ends at 105°, making its pulse width 30° (their 
difference); its three regions of varying peak-to-peak 
voltage (and, hence, differing %V pplavg ) are the intervals of 
a from 0° to 15°, 15° to 75°, and 75° to 105°. 

Table 3. Angular characteristics of some p-pulse 

controlled rectifiers with purely resistive loads. 



0L 

6u 

pulse-width 

® limit 1 

Otlimit 2 

^upper limit 

pulse-count 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

6 

60 

120 

60 

30 

60 

120 

12 

75 

105 

30 

15 

75 

105 

18 

80 

100 

20 

10 

80 

100 

24 

82.5 

97.5 

15 

7.5 

82.5 

97.5 

96 

88.125 

91.875 

3.75 

1.875 

88.125 

91.875 


Just as Table 1 was developed to provide L crit for 
selected values of a , Table 4 does the same for %V pp/avg for 

the waveforms of Figure 19. Firing angles of 120° for the 
6-pulse controlled rectifier and 100° for the 18-pulse 
controlled rectifier would result in infinite %V /avg because 

those values of a are the upper limits (and, therefore, not 
inclusive) for the two rectifiers, in accordance with Table 

3. 
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Table 4. 


Peak-to-peak ripple voltage as a percentage of 
the average output voltage for selected firing 
angles for the controlled p-pulse rectifier 
waveforms of Figure 19. 



firing angle alp 

la (degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

6 

14% 

25% 

40% 

60% 

88% 

125% 

181% 

225% 

288% 

391% 

594% 

1197% 

12 

3% 

10% 

19% 

30% 

44% 

62% 

91% 

144% 

236% 

398% 

1199% 


18 

2% 

6% 

13% 

20% 

29% 

42% 

60% 

96% 

198% 

399% 

OO 


24 

<1% 

5% 

10% 

15% 

22% 

31% 

45% 

72% 

148% 

399% 



96 

<1% 

1% 

2% 

4% 

5% 

8% 

11% 

18% 

37% 

400% 




While the data presented by Table 4 is somewhat useful 
for quickly finding the value of %V pplavg for a given value of 

a for a particular p-pulse controlled rectifier with a 
resistive load, another way to look at the data is by 
assigning the %V pplavg of the 6-pulse controlled rectifier as 

the basis to which the others are compared. The result is 
Table 5. For realistic values of a, higher pulse count 
controlled rectifiers clearly possess much less %V pp/avg than 
lower pulse count variants. 
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Table 5. 


Peak-to-peak ripple voltage as a percentage of 
the average output voltage as a percentage of 6- 
pulse %V /av for selected firing angles for the 

controlled p-pulse rectifier waveforms of Figure 

19. 




iring angle a 

pha (c 

egrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

6 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

25% 

39% 

48% 

50% 

50% 

50% 

50% 

64% 

82% 

102% 

202% 


18 

11% 

25% 

32% 

33% 

33% 

33% 

33% 

43% 

69% 

102% 

oo 


24 

6% 

19% 

24% 

25% 

25% 

25% 

25% 

32% 

52% 

102% 



96 

<1% 

5% 

6% 

6% 

6% 

6% 

6% 

8% 

13% 

102% 




B. OUTPUT VOLTAGE RIPPLE FOR A CONTROLLED RECTIFIER WITH 
A TWO-POLE LC OUTPUT FILTER 


Turning back to the circuit of Figure 9, we now desire 
to calculate the output ripple for the type of circuit we 
are primarily interested in analyzing. It turns out that 
the voltage ripple calculations for a controlled rectifier 
with a two-pole LC output filter are actually somewhat 
simpler than those for the controlled rectifier with a 
purely resistive load because we are going to assume the 
inductor is always large enough to maintain a continuous 
current; therefore, there is only one case to evaluate, not 
two. Additionally, some of the calculations are the same as 
for the purely resistive case, so we need not repeat our 
previous work. To begin from [4], the average capacitor 
voltage is given by (70), and the RMS load voltage is given 
by (134) . Then, the RMS ripple voltage is given by 




RMS-ripple 




2 

C(avg) 


(153) 
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Here we will make a small notation change, letting V DC be 
equal to V avg and V AC be equal to V RMS _ rippIe to better illustrate 
the relationship between RMS, DC, and AC quantities. 
Additionally, it is important to understand that V AC is 
applied across the entire load (which includes the 
resistor, capacitor, and the inductor). Therefore, (153) 
becomes 

y AC =4Vms-Vi c (154) 

and (141) becomes 


v = V — 

v AC v m\\ 2 


1 + 


Fpcos(2a ) 


2 n 


A 2 p 2 2 

-— cos a 

n 


'155; 


The impedance of the capacitor is 

Z c s —-— 
jpcoC 


(156) 


and the impedance of the inductor is 

Z L = jpcoL , (157) 

where these are considered approximations because Z c and Z L 
are composed of the frequencies pco , 2 p(D, 3pco , etc.; 

however, the dominant frequency is pco . Next, the impedance 
of the capacitor in parallel with the resistor is 


7 - 

^R\\C 


R 


1+ jpcoRC 


(158) 


which can be written in terms of real and imaginary parts 
as 


7 - 

^R\\C ~ 


R 


pcoRC 


1 + (pcoRCy 1 + (pcoRC) 


(159) 
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The impedance of Z L in series with Z fi||c is 


7 - 

total ~ 


R 


\ +(pcoRC)' 

which can be rewritten as 




r pcoRC 
pcoL - T 

1 + {pooRC) 


7 - 

total ~ 


R 


+ ./' 


pcoL + p 3 oo'RLC 2 - pooR 2 C 


l + (pcoRC ) 

The magnitude of Z tofa/ becomes 


l + (/^C) 2 


7 - 

total ~ . 


r r 

2 

_l_ 

pcoL + p 'oj R LC^ - pcoR 1 C 

[1 + ( P coRC) 2 

1 

1 + (pcoRC ) 2 


We can now solve for / AC from Ohm's Law as 


i , V 

I J ~ V AC 

KACl— y 

^total 


By substituting (155) and (162) into (163), we get 


I AC = 



r r 

2 

pcoL + p 3 co 3 R 2 LC 2 - pcoR 2 C 

[1 + (pcoRC) 2 

1 

1 + (pcoRC) 2 


Therefore, the capacitor RMS voltage V C{RMS) is given by 


V ~ J 7 

v C(RMS) ~ 1 AC C 


and the capacitor peak-to-peak ripple voltage is 


^C(pp-ripple) — (2V2 jV c(SMS) 


(160) 


(161) 


(162) 


(163) 


(164) 


(165) 


(16 6 ) 
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Here, it is important to note another approximation that 


has taken place. We have made the assumption that V C{RMS) is a 
perfect sinusoid such that its peak value is simply equal 
yf2\f 

to c(rms) _ Finally, the capacitor peak-to-peak ripple 

voltage as a percentage of the average output voltage 

0 // ° V C(pp/ a vg) 1 s 


%V, 


C(pp/avg ) 


C(pp-ripple) 

V "c(avg) J 


100 % , 


(167) 


where V C(avg) is given by (70) and V c(pp _ ripple) is given by (166) . 

One of the things that would be interesting to 
investigate is how %V„, , s behaves as the individual 
circuit elements (inductor, capacitor, and resistor) are 
changed one at a time. While the behavior of 0 ^ C (ppiavg) f° r 
every circuit is unique depending on its peak line-to-line 
input voltage V m , pulse count, and the values of the 
inductor, capacitor, and resistor, we will examine one case 
with the aid of Figure 20. The particular circuit that 
generated the waveforms of Figure 20 is a 6-pulse 
controlled rectifier (in the configuration given by Figure 
9) with a V m of 100 V, a 25 mH inductor, a 50 capacitor, 

and a 10Q resistor. Through some initial trial and error 
with different resistance values, it was discovered that 
°/oV C (ppiavg) was ver Y insensitive to changes in the resistance, 

so it was excluded from further sensitivity analysis. The 
initial values of the inductor and capacitor are as 
previously mentioned and are designated as L 0 and C 0 in the 
legend of Figure 20. 
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Figure 20. Capacitor peak-to-peak ripple voltage 

as a percentage of the average output voltage for a 6- 
pulse controlled rectifier with a two-pole LC output 

filter. 

Somewhat arbitrarily, it was decided to increase and 
decrease L 0 and C 0 by fifty percent to observe the effect on 
%V r , , , (care had to be taken to ensure the circuit was 

C( pp/avg J 

still in CCM for the reduced inductor value) . The larger 
inductor and capacitor are designated as L H and C H in the 
legend of Figure 20, while the reduced components are L l 
and C L , respectively. Interestingly, the combination of a 
reduced inductor and initial capacitor (worst %^qpp/ ors )) 

produced close results to the initial inductor and reduced 
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capacitor pair, while the larger inductor and initial 
capacitor produced nearly identical results to the 
combination of the initial inductor and bigger capacitor— 
both of which could be said to have the best °/oV C ( pp/avg ) • 

In order to see how increasing pulse count affects 
°/oV C (p p/avg y 0// °Vc(p P /avg) was plotted in Figure 21 for a 6-, 12-, 

18-, 24-, and 96-pulse controlled rectifier with a two-pole 

LC output filter using the same initial component values 
and V„ as before. 

m 



a (degrees) 


Figure 21. Capacitor peak-to-peak ripple voltage 

as a percentage of the average output voltage for a 
6-, 12-, 18-, 24-, and 96-pulse controlled rectifier 
with a two-pole LC output filter. 


60 

























Much like the plots of Figures 15 and 19, tremendous 
improvement is noted in increasing from a 6-pulse to a 12- 
pulse controlled circuit, but further increasing pulse 
count results in less of an improvement. 

The values of %V c ^ pp/avg ^ for selected values of a for the 

waveforms of Figure 21 are collected in tabular form in 
Table 6. The capacitor peak-to-peak ripple voltage as a 
percentage of the average output voltage approaches 
infinity at the firing angle limit of 90° because } is 

zero at 90°; hence, the firing angle limit is exclusive of 
90° . 


Table 6. Capacitor peak-to-peak ripple voltage as a 
percentage of the average output voltage for 
selected firing angles for the controlled p-pulse 
rectifier waveforms of Figure 21. 



firing angle alpha (degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

2% 

3% 

6% 

9% 

13% 

18% 

26% 

41% 

84% 

12 

<1% 

<1% 

<1% 

1% 

1% 

2% 

3% 

5% 

10% 

18 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

1% 

3% 

24 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

1% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 


As stated previously, every circuit of this type is 
unique and will produce differing %V c ^ /av ^ depending on the 

parameters of the controlled rectifier with a two-pole LC 
output filter; therefore. Table 6 has limited utility. 
However, if we employ the technique we have previously used 
and treat 0 ^ C (p P iavg) f° r the 6-pulse controlled rectifier as 
the basis to which the other rectifiers are compared, we 
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generate Table 7. As opposed to Table 6, which deals in 
absolutes, the relationships in Table 7 are always valid no 
matter what the controlled rectifier circuit parameters and 


value of V m are (so long as the circuit remains in CCM). 


Table 7. Capacitor peak-to-peak ripple voltage as a 
percentage of the average output voltage as a 
percentage of 6-pulse %V c ^ /av ^ for selected firing 

angles for the controlled p-pulse rectifier 
waveforms of Figure 21. 



firing angl 

le alpha (degrees) j 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

5.8% 

9.9% 

11.1% 

11.4% 

11.6% 

11.6% 

11.6% 

11.7% 

11.7% 

18 

1.1% 

2.8% 

3.2% 

3.3% 

3.4% 

3.4% 

3.4% 

3.4% 

3.4% 

24 

<1% 

1.1% 

1.3% 

1.4% 

1.4% 

1.4% 

1.4% 

1.4% 

1.4% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 
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V. 


HARMONICS 


A. BACKGROUND 

Harmonics have existed on power systems from the time 
of the very first generators. At that time, the harmonic 
components were so small that their effects on systems were 
negligible. This was due to the lack of non-linear loads 
before the 1960s [12] . However, in recent years, there has 
been a significant increase in the installation and use of 
power electronics equipment onboard ships and at offshore 
installations. The operation of this equipment can, and has 
in many cases, significantly degraded the ship or offshore 
installation electrical power quality to such an extent 
that corrective measures have to be implemented in order to 
reduce the resultant adverse effects on the electrical 
plant and equipment. The quality and security of voltage 
supplies are important to the safety of any vessel and its 
crew and to the protection of the marine environment. Any 
failure or malfunction of equipment such as propulsion or 
navigation systems can result in an accident at sea or 
close ashore with dire consequences. What we are talking 
about is harmonic distortion of voltage supplies caused by 
the operation of electronic devices which draw non-linear 
(i.e., non-sinusoidal in nature) currents from the voltage 
supplies. These very same pieces of non-linear equipment 
can also be adversely affected by harmonic currents and the 
subsequent voltage distortion they produce, as can the 
majority of linear equipment (particularly generators, AC 
motors and transformers) [13]. 
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B. COMMON PROBLEMS ASSOCIATED WITH HARMONICS 

Power converters are considered the primary source of 
undesired harmonics. These AC/DC power converters are 
extensively used in various applications like power 
supplies, DC motor drives, front-end converters in 
adjustable-speed AC drives, high voltage DC transmission, 
switched-mode power systems, fluorescent lights (with 
electronic ballasts), utility interfaces with non- 
conventional energy sources, in-process technology like 
welding, power supplies for telecommunications systems, 
aerospace, and military environments. Traditionally, AC/DC 
power conversion has been dominated by diode or phase 
controlled rectifiers which act as non-linear loads on 
power systems and draw input currents that are rich in 
harmonics and have poor supply power factor, thereby 
creating the power quality problem for the distribution 
network and for other electrical systems in the same 
vicinity of the rectifier. Some common problems associated 
with the production of these harmonics include: 

• Large reactive power drawn by rectifiers from the 
power system which requires that the distribution 
equipment handle large power, thus, increasing 
its volt-ampere ratings; 

• Voltage drops at the buses; 

• Higher input current harmonics resulting in the 
distorted line current which tends to distort the 
line voltage waveform. This often creates 
problems in the reliable operation of sensitive 
equipment operating on the same bus; 

• Increased losses in the equipment such as 
transformers and motors connected to the utility; 

• Electromagnetic interference with any nearby 
communications circuits; 
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• Excessive neutral current, resulting in 
overheated neutrals. The currents of triplen 
harmonics, especially odd harmonics (3 rd , 9th 
15th, etc.) are actually additive in the neutral 
of three-phase, wye-connected circuits; 

• Measuring equipment and meters reading 
incorrectly, including induction disc-type Watt- 
hour meters and averaging type current meters; 

• Blown fuses on power factor correction capacitors 
due to high voltages and currents from resonance 
with line impedance and capacitor bank failures; 

• Mal-operation of equipment such as computers, 
telephone systems, and electronic controllers; 

• Nuisance operation of protective devices 
including false tripping of relays and failure of 
uninterruptible power supplies to transfer 
properly, especially if the controls incorporate 
zero-crossing sensing circuits; and 

• Damaging dielectric heating in cables [14]. 

C. VISUAL REPRESENTATIONS OF HARMONICS 

Before delving directly into the mathematics and 
equations dealing with harmonics, it is useful to go a bit 
deeper into what is meant by a load being either linear or 
non-linear. A linear load will draw current that is 
proportional to the applied voltage, as in Figure 22 [13]. 



Voltage and current waveforms for a 
linear load. From [13]. 
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Figure 22. 





It stands to reason then, that a non-linear load will draw 
a current that does not flow in proportion to the applied 
voltage [15] . One such example of a non-linear load is the 
load current profile for a three-phase AC drive shown in 
Figure 23. 



Figure 23. Input line to neutral voltage (in black) 

and input current (in green) of phase A for a three- 
phase AC drive. From [15]. 

Around the 1830s, a French mathematician named Fourier 
discovered that a distorted waveform—like the current 
waveform shown in Figure 23—can be represented as a series 
of sine waves, each an integer number multiple of the 
fundamental frequency and each with a specific magnitude. 
For example, for a 60 Hz system, the fundamental frequency 
is 60 Hz and the 5th harmonic would be five times 60 Hz or 
300 Hz [15] . When all harmonics are added to the 
fundamental, a waveform known as a complex wave is formed. 
A simple example of a complex wave consisting of the 
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fundamental and 3 rd harmonic is shown in Figure 24. The 
complex wave is clearly the algebraic summation of the 
other two waves [13] . 

60 Hz Fundamental 



Figure 24. Construction of a complex wave from 

the fundamental and 3 rd harmonic. From [13]. 

The collective sum of the fundamental and each 
harmonic is called a Fourier series. This series is often 
viewed as a spectrum analysis where the fundamental 
frequency and each harmonic component are displayed 
graphically in a bar chart format as shown in Figure 25. 
The total current is the square root of the sum of the 
squares of the individual harmonic components. One general 
thing we can say regarding the harmonic content of 
something is by simply looking at the wave shape. The 
more it looks like a sine wave, the lower the harmonic 
content [15]. 
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Harmonic components 

0.2 
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Fundamental 5th 7th 11th 13th 17th 19th 


Figure 25. Harmonic spectrum analysis. From [15]. 


D. FOURIER SERIES AND TOTAL HARMONIC DISTORTION 

From [16], in general, any non-sinusoidal waveform f(t ) 
repeating with an angular frequency Q) can be expressed as 

OO 00 

/«) = F„ + X/,«) = t<i„+Z (a h cos (hcot) + b h sin(/i<y?)}, (168) 

h =1 ^ h =1 


where h is the harmonic number and F 0 = a 0 /2 is the average 
value. Additionally, we have 


and 


1 2 f 

a h =—\ f (t) cos (hcot)d[cot) 
n i 


(169) 


K 


1 2 r 

— f (t) sin (hcot^d(at) , 

n o 


(170) 


in which h = 0,...,oo . From (168) and (169), we can solve for 
the average value F 0 (taking note that a = 2n!T) as 


\_ 

2 


= —a n = J f(t)d (at) = ^j[ f(t)dt 


2 n 


T 


(171) 


Oftentimes, AC waveforms have a zero average value. 
Moreover, by the use of waveform symmetry, it is often 
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possible to simplify the calculations of a h and b h in (169) 
and (170) . A list of the types of symmetry, requirements, 
and formulas for a h and b h is provided in Table 8. 


Table 8. Use of symmetry in Fourier analysis. From [16]. 


Symmetry 


Condition Required 


a h and b h 


Even 

Odd 

Half-wave 


Even 

quarter-wave 


/(—o =m 


fi-t ) = -At) 


fit) = -fit + \T) 


Even and half-wave 


Odd Odd and half-wave 

quarter-wave 


2 fir 

bh = 0 a h = - f(t)cos(hoit) d(ujt) 

it Jo 

2 fir 

Oh = 0 b h = — /(r)sin(Ao)/) d(ut) 

Tt Jo 

a h = b h - 0 for even h 

2 fit 

a h = — f f(t) cos(hoit) d(<ot) for odd h 
Jo 

2 r tt 

fe/, = — /(t)sin(/uo0 d((or) for odd h 
tt Jo 


b h = 0 for all h 


a h 


= 


4 fir/2 

TT J 0 " 


f {t)cos(hoit) d(oit) for odd h 
for even h 


a h = 0 for all h 


I 4 fir/2 

_ I — J f(t)sin(hoit) d(oit) for odd h 


for even h 


For simplicity, let us now assume that we have a 
purely sinusoidal utility input voltage at the fundamental 

frequency (with co x =(0 and f\=f) as v s (f) = yp2V s sin&J/ . Then, 

from [16], the input current in steady state is the sum of 
its Fourier (harmonic) components, given by 

U0 = 4i(0 + Z^(0' < 172 ) 

h^\ 
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where we have assumed there is no DC component in i s , i sl is 
the fundamental (line frequency ) component, and i sh is the 
component at the h harmonic frequency f h ( = hfi)- T ^e amount of 

distortion in the current waveform is quantified by means 
of an index termed the THD. The distortion part of the 
current i dis is given by 

0) = K (0 - 4i (0 = X l sh (0 • ( 17 3 ) 

h* 1 

In terms of RMS values, we have 

' (!74) 

where the RMS value of i s (t) has been determined by 


! s * T 


M[4(0] 2 ^ 


(175) 


Finally, the THD of the current is defined as 


%77/D, = 100% 


r I ^ 

1 du 

\ J sl J 


= 100 % 
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E. REDUCTION OF HARMONICS BY INCREASING PULSE NUMBER 

One of the most important benefits gained by 
increasing the pulse count of a controlled rectifier is a 
reduction in harmonic content. Without diving into the 
detailed mathematics about why this is the case, suffice it 
to say that when multiple properly phase shifted converters 
are connected in series or parallel, harmonic cancellation 
occurs such that the line currents drawn by a p-pulse 
converter system have a harmonic content given by 
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h = kp± 1 


(177) 


where k = 1,2,3... [6]. Additionally, from [16], the RMS values 
of the harmonic components are inversely proportional to 
their harmonic order as 



From (177) and (178), it is easy to understand that higher 
pulse count systems will have smaller harmonics that are 
further away from the fundamental component. 

It is natural to consider a 6-pulse controlled 
rectifier system as the basic building block from which 
higher pulse count systems are constructed. In order to 
achieve the desired harmonic elimination using stacked 6- 
pulse rectifiers, the required phase shift between 
individual rectifiers is given as phase shift (f> equals 60° 
divided by the number of 6-pulse converters [16]. For 
example, a 12-pulse converter would be made from two 6- 
pulse converters, and 60°/2=30°; therefore, the required 
phase shift between the two 6-pulse converters is 30°. One 
realization of such a system is shown in Figure 26. The 
system consists of two bridge rectifiers connected so that 
the DC outputs are in series and the two voltages add. The 
AC inputs to the two bridges are isolated by two 
transformers that are connected wye-delta and delta-delta, 
respectively. The transformers are connected to take 
advantage of the differing phase shifts of the two 
connections [3]. 
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Figure 26. Series circuit arrangement for a 12-pulse 

rectifier using a wye-delta and delta-delta connected 
transformer to achieve the necessary 30° phase shift 
between 6-pulse rectifier units. From [3]. 

While it was relatively simple to realize the 12-pulse 
rectifier by taking advantage of the 30° phase shift that 
naturally occurs for a wye-delta transformer, how would we 
go about achieving other phase shifts? For example, for an 
18- and 24-pulse count rectifier, we would require phase 
shifts of 20° and 15°, respectively. One such device that 
can easily accomplish this is called the zig-zag 


transformer [17] . 




VI. SIMULINK® SIMULATIONS 


In order to validate the relationship between pulse 
count and harmonic content (including THD), simulations for 
a 6-, 12-, 18-, and 24-pulse controlled rectifier with a 

two-pole LC output filter were conducted using Simulink®, 
specifically looking at the harmonic content of the phase A 
input line current i . Additionally, it was desired to 

observe the behavior of the harmonics and %THD. as a was 
varied. 

A. CIRCUIT PARAMETERS COMMON TO ALL SIMULATION MODELS 

The fixed step-size for all simulations was 10 |a.s 
using the Runge-Kutta (ode4) solver; total simulation time 
was 1.0 s. 

The voltage source was a three-phase 60 Hz ramp source 
with a 2, 980 V/s slope, going from zero to effectively 150 
V (RMS, line-to-line) in 0.05 s. The use of a ramp source 
rather than a step source allows for faster steady-state 
results. This voltage source was chosen to match that of an 
example problem from [4]. 

The transformers were zigzag phase shifting 
transformers with the following common parameters: 

• Secondary winding (abc) connection: Y, 

• Nominal power: 1 MVA, 

• Nominal frequency: 60 Hz, 

• Primary nominal voltage V p : 150 V (RMS, line-to- 

line) , 

• Winding 1 resistance : 0, 
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• Winding 1 inductance L, : 0, 

• Winding 2 resistance R 2 : 0, 

• Winding 2 inductance L 2 : 0, 

• Winding 3 resistance R 3 : 0, 

• Winding 3 inductance L 3 : 0, 

• Magnetizing branch resistance R m : 10 kf2 , and 

• Magnetizing branch inductance L m : "inf." 

The main circuit component values were as follows: 

• Inductor: 23.375 mH, 

• Capacitor: 200 hF, and 

• Resistor: 19.3 Q. 

These values were also chosen to match those of the same 
example problem from [4]. 

For the discrete synchronized 6-pulse generators, a 
pulse width of 10° was used with double-pulsing disabled. 

For the 6-pulse controlled rectifier module, an 
inductance of 1 pH was added to each branch feeding into 
the universal bridges. 

The universal bridges had the following common 
parameters: 

• Number of bridge arms: 3, 

• Snubber resistance R s : 100 Q, 

• Snubber capacitance C s : 0.1 |iF, 

• Power electronic device: thyristors, 

• On-resistance R on : 1 mQ , 

• On-inductance L on : 0, and 

• Forward voltage V f : 1 mV. 
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The Simulink® schematic for the internal view of the 
6-pulse controlled rectifier module common to all of the p- 
pulse simulation models is shown in Figure 27. 



Vc 


Figure 27. Simulink® schematic of internal view of 
6-pulse controlled rectifier module. 

B. DATA COLLECTION 

As stated previously, the data of interest was the 
harmonic content of i as a varied. Because it was 
desirable to observe how individual harmonics (normalized 
to the fundamental) and %THD i varied with respect to a, it 
was decided to evaluate the first ten non-zero harmonics 
(not including the fundamental) for each p-pulse rectifier. 
The first ten non-zero harmonics for the simulated p-pulse 
rectifiers are shown in Table 9. 
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Table 9. Harmonic numbers for the first ten non-zero 

harmonics (not including the fundamental) for some 

p-pulse rectifiers. 


6-pulse 

12-pulse 

18-pulse 

24-pulse 

5 

11 

17 

23 

7 

13 

19 

25 

11 

23 

35 

47 

13 

25 

37 

49 

17 

35 

53 

71 

19 

37 

55 

73 

23 

47 

71 

95 

25 

49 

73 

97 

29 

59 

89 

119 

31 

61 

91 

121 


In order to obtain the RMS values of the individual 
harmonics, a string of eleven Fourier analyser blocks was 

used in series with gain blocks set at V2/2, specifying the 
fundamental frequency (60 Hz) and harmonic number of 
interest for each Fourier analyser block. Additionally, to 
obtain %THD. , a total harmonic distortion block in series 
with a gain block equal to 100 was used, again specifying 
the fundamental frequency of 60 Hz. The Simulink® schematic 
symbology depicting this is shown in Figure 28. 
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Figure 28. Obtaining the RMS values for the 
fundamental and first ten non-zero harmonics and %THD i 

from Simulink®. 


Because it was initially thought that 
and %THD i would change faster at lower 
(erroneously, in this particular case), the 
size for a was chosen to be 3°, starting at 


the harmonics 
values of a 
initial Step¬ 
an 0° and going 


77 










































































































then increasing step-size to 10 


and 


up to a = 30 °, 
continuing up to a = 80 °. 
a = 90 ° (or even a = 80 ° for 
L crit would approach inf ini 
would no longer be in CCM. 


It was not possible to go up to 
all controlled rectifiers) since 
ty, and the controlled rectifier 


The ten harmonics for each a were then normalized to 
the fundamental by dividing the RMS value of each harmonic 
by the RMS value of the fundamental. The normalized 
harmonics were then plotted versus a for each p-pulse 
controlled rectifier. Additionally, the %THD j for all four 
of the p-pulse controlled rectifiers was plotted versus a . 


C. CONTROLLED 6-PULSE RECTIFIER SIMULATION 

The controlled 6-pulse rectifier Simulink® model is 
shown in Figure 29. 



Figure 29. Simulink® model for a controlled 6-pulse 
rectifier with a two-pole LC output filter. 

Although a transformer is not necessary for the controlled 
6-pulse rectifier (because no other phase shifts beyond 
those that already exist for the three-phase source are 
required), a zigzag phase shifting transformer was used to 
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provide isolation. The secondary nominal voltage V 3 and 
phase shift (j) were 150 V (RMS, line-to-line) and 1X10- 
degrees (because 0° is not permitted), respectively. 

The first ten normalized harmonics are plotted versus 
a in Figure 30. The 6-pulse controlled rectifier did not 
have enough inductance at a = 80 ° to remain in CCM. 
Interestingly, the first normalized non-zero harmonic ( h 5 ) 

increases with a (which causes %THD j to worsen) , but the 
second non-zero harmonic {h 7 ) decreases with a; the 
remaining eight normalized harmonics appear to only 
decrease very slightly as a increases. 



Figure 30. First ten non-zero normalized harmonics 

versus a for a 6-pulse controlled rectifier with a 
two-pole LC output filter. 
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D. CONTROLLED 12-PULSE RECTIFIER SIMULATION 

The controlled 12-pulse rectifier Simulink® model is 
shown in Figure 31. 



6-pulse controlled 
rectifier 2 

Figure 31. Simulink® model for a controlled 12-pulse 
rectifier with a two-pole LC output filter. 

The secondary nominal voltage V 3 for both the upper and 
lower zigzag phase shifting transformer were 75 V (RMS, 
line-to-line) while the phase shift $ was 15° for the upper 
transformer and -15° for the lower transformer. The reason 
the secondary nominal voltages are the same and are half of 
what they were for the 6-pulse case is because it is 
desirable to maintain the same overall rectifier output 
side voltage. Recall that since the two 6-pulse controlled 
rectifiers are connected in series, their output voltages 
are additive. The ± 15 ° phase shifts create the necessary 30° 
phase difference between the two 6-pulse controlled 
rectifiers to ensure proper harmonic cancellation. 
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The first ten normalized harmonics are plotted versus 
a in Figure 32. The normalized harmonics of the 12-pulse 
controlled rectifier displayed fairly similar behavior to 
those of the 6-pulse controlled rectifier, although, they 
are smaller, as predicted by (177) and (178) . 



0 10 20 30 40 50 60 70 80 


a (degrees) 

—♦— hll -B-hl3 -±-h23 -*-h25 -*-h35 
-#-h37 h47 h49 h59— ♦— h61 


Figure 32. First ten non-zero normalized harmonics 

versus a for a 12-pulse controlled rectifier with a 
two-pole LC output filter. 


E. CONTROLLED 18-PULSE RECTIFIER SIMULATION 

The controlled 18-pulse rectifier Simulink® model is 
shown in Figure 33. 
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6-pulse controlled 
rectifier 3 


Figure 33. Simulink® model for a controlled 18-pulse 
rectifier with a two-pole LC output filter. 

The secondary nominal voltage V 3 for all three zigzag 
phase shifting transformer was 50 V (RMS, line-to-line) 
while the phase shift ^ was 20° for the upper transformer, 
1X10- ; degrees (because 0° is not permitted) for the middle 
transformer, and -20° for the lower transformer. The reason 
the secondary nominal voltages are the same is for the same 
reason as before. The ± 20 ° and 0° (effectively) phase shifts 
create the necessary 20° phase differences between the 
three 6-pulse controlled rectifiers to ensure proper 
harmonic cancellation. 

The first ten normalized harmonics are plotted versus 
a in Figure 34. The normalized harmonics of the 18-pulse 
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controlled rectifier displayed similar behavior to those of 
the other two controlled rectifiers, although, they are 
smaller, again as predicted. 



Figure 34. First ten non-zero normalized harmonics 

versus a for an 18-pulse controlled rectifier with a 
two-pole LC output filter. 


F. CONTROLLED 24-PULSE RECTIFIER SIMULATION 

The controlled 24-pulse rectifier Simulink® model is 
shown in Figure 35. 
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6-pike coitrolled 


Figure 35. Simulink® model for a controlled 24-pulse 
rectifier with a two-pole LC output filter. 

The secondary nominal voltage V 3 for all four zigzag 
phase shifting transformer was 37.5 V (RMS, line-to-line) 
while the phase shift ^ was 22.5° for the top transformer, 
7.5° for the second transformer, -7.5° for the third 
transformer, and -22.5° for the bottom transformer. The 
phase shifts of + 22 . 5 ° and ± 7 . 5 ° create the necessary 15° 
phase differences between the four 6-pulse controlled 
rectifiers to ensure proper harmonic cancellation. 
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The first ten normalized harmonics are plotted versus 
a in Figure 36. It is no surprise that the normalized 
harmonics of the 24-pulse controlled rectifier displayed 
similar behavior to those of the other three controlled 
rectifiers, although, they are once again smaller, just as 
predicted. 



—h23 
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_±_h47 -*-h49 
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-*-h73 
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Figure 36. First ten non-zero normalized harmonics 

versus a for a 24-pulse controlled rectifier with a 
two-pole LC output filter. 


G. SIMULATION RESULTS FOR %THD i 

The %THD j versus a for p-pulse controlled rectifiers 
with a two-pole LC output filter is shown in Figure 37. 
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Figure 37. Percent total harmonic distortion in i 

versus a for p-pulse controlled rectifiers with a two- 

pole LC output filter. 


The data clearly show that a higher pulse count produces a 
lower %77/D, , with the 24-pulse controlled rectifier the 
clear winner. However, an incorrect conclusion that might 
be drawn from Figure 37 is that %THD i does not vary much 
with a ; in fact, it appears as if %THD i only gets slightly 
worse for very high values of a . The real reason that 
%THD i does not appreciably worsen as a increases is 
because the value of the inductor chosen for the 
simulations was very large (23.375 mH) . The point is that 
%THD i and L crit are both related to a in much the same way. 

, Figure 37 is related to Figure 15. If we had 
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instead initially selected a much smaller inductor that 
would still ensure CCM for the 6-pulse controlled rectifier 
up to a- 30° ( = 3 mH) , we would more than likely see %THD i 

increase for much smaller values of a , particularly for the 
6-pulse controlled rectifier and less so for the other 
higher pulse-count controlled rectifiers because of the 
reduction in L crjt as pulse count increases. 
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VII. CONCLUSIONS 


A. SUMMARY OF RESULTS 

1. Critical Inductance 

The critical inductance as a function of a was 
calculated for a p-pulse controlled rectifier with a two- 
pole LC output filter. Plotting L crjt for several multi-pulse 
controlled rectifiers demonstrated how it increased as a 
increased and that higher pulse count controlled rectifiers 
have a lower L crit than lower pulse count controlled 
rectifiers for a given value of a . Additionally, L crit as a 
percentage of L crit for a 6-pulse controlled rectifier was 
calculated for selected values of a (0° to 90°, in 10° 
increments) for a 6-, 12-, 18-, 24-, and 96-pulse 

controlled rectifier. This clearly illustrated how the 
greatest reduction in L crit occurs in increasing from a 6- 
pulse to a 12-pulse controlled rectifier and that, while 
further increasing pulse count does reduce L crit , the 
reduction is less. 

2. DC Output Voltage Ripple 

The peak-to-peak ripple voltage as a percentage of the 
average output voltage as a function of a was calculated 
for a p-pulse controlled rectifier with a purely resistive 
load. Plotting %V pp/avg for several multi-pulse controlled 

rectifiers demonstrated how it increased as a increased and 
that higher pulse count controlled rectifiers have a lower 
%Vp p , aV g than lower pulse count controlled rectifiers for a 
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given value of a . Additionally, %V f as a percentage of 
%Vp p / m , g for a 6-pulse controlled rectifier was calculated for 

selected values of a (0° to 90°+, in 10° increments) for a 
6-, 12-, 18-, 24-, and 96-pulse controlled rectifier. This 

clearly illustrated how the greatest reduction in %V pplavg 

occurs in increasing from a 6-pulse to a 12-pulse 
controlled rectifier and that, while further increasing 
pulse count does reduce %V , , the reduction is less. 

The peak-to-peak ripple voltage as a percentage of the 
average output voltage as a function of a was also 
calculated for a p-pulse controlled rectifier with a two- 
pole LC output filter. Plotting %V pp/avg for several multi¬ 
pulse controlled rectifiers demonstrated how it increased 
as a increased and that higher pulse count controlled 
rectifiers have a lower %V pplavg than lower pulse count 
controlled rectifiers for a given value of a . Additionally, 
%V pp ,avg as a percentage of %V pplavg for a 6-pulse controlled 

rectifier was calculated for selected values of a (0° to 
90°, in 10° increments) for a 6-, 12-, 18-, 24-, and 96- 

pulse controlled rectifier. This clearly illustrated how 
the greatest reduction in %V pplavg occurs in increasing from a 
6-pulse to a 12-pulse controlled rectifier and that, while 
further increasing pulse count does reduce °/oV pp/avg , the 
reduction is less. 

3. Harmonics and Total Harmonic Distortion 

Simulations were conducted using Simulink® for a 6-, 

12-, 18-, and 24-pulse controlled rectifier with a two-pole 

LC output filter, specifically investigating the behavior 
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of the first ten non-zero harmonics (normalized to the 
fundamental) and %THD j of i for specific values of a. 

Plotting the harmonics of each p-pulse controlled rectifier 
validated that the harmonics for higher pulse count 
controlled rectifiers are smaller than those for lower 
pulse count controlled rectifiers. Additionally, it 
revealed that the first harmonic increased as a increased, 
while the second harmonic decreased as a increased. The 
remaining eight harmonics did not show significant change 
as a varied. 

Plotting %THD i for the p-pulse controlled rectifiers 
confirmed that it increased as a increased (though it was 
reasoned that it would show a more definitive increase as a 
increased if a smaller inductor had been used) and that 
higher pulse count controlled rectifiers have a lower %THD i 
than lower pulse count controlled rectifiers for a given 
value of a . This also illustrated how the greatest 
reduction in %THD i occurs in increasing from a 6-pulse to a 
12-pulse controlled rectifier and that, while further 
increasing pulse count does reduce %THD j , it results in 
diminishing returns. 

B. RECOMMENDATIONS FOR FUTURE WORK 

1. Discontinuous Conduction Mode 

When the inductor in a multi-pulse phase controlled 
rectifier with a two-pole LC output filter is less than the 
value of L crit , the inductor current is no longer continuous, 
and the equations, tables, and figures of Chapters III and 
IV no longer apply. Although operation of a controlled 
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rectifier in discontinuous conduction mode is not generally 
considered desirable, it is still a possible mode for the 
rectifier, and the model is not complete without it. 

2. Size and Weight Considerations 

While it is clear that a higher pulse count controlled 
rectifier produces better results than one with a lower 
pulse count, it can only do so at the cost of greater 
system complexity (and hence, less system reliability) and 
with more power electronic components. Power SCRs, 
inductors, capacitors, and especially low frequency 
transformers are bulky and heavy; this must be considered 
in real systems designed to operate on Navy ships with 
space and weight limitations. One useful study would be to 
use actual space and weight characteristics of the power 
electronic components for these controlled rectifiers and 
to determine the optimal pulse count system that would 
satisfy the electrical requirements and be within size and 
weight constraints for the intended Navy ship platform. 

3. Parallel Vice Series Connection of Multi-Pulse 
Controlled Rectifiers 

Only the series (vice parallel) connection of multiple 

controlled rectifiers in order to increase pulse count was 

considered in this thesis. With a series connection, the 

output voltage of the system is the summation of the 

individual output voltages of each rectifier, whereas, with 

a parallel connection, the output current of the system is 

the summation of the individual output currents of each 

rectifier. However, connecting rectifiers in parallel is 

more complex due to fluctuations in DC voltages which cause 

circulating currents to flow. To remedy this, these types 
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of systems typically require the use of an interphase 
transformer which absorbs the differences in the DC 
voltages between individual rectifiers. A useful study 
would be the evaluation of differences between the two 
topologies as a result of the added interphase transformer. 

4. The Effect of a on Total Harmonic Distortion 


In general, increasing a 

will 

cause 

the THD of 

a 

controlled 

rectifier system 

to 

worsen. 

However, 

a 

theoretical 

mathematical model 

for 

how THD 

behaves as 

a 


function of a similar to how L crit behaves as a changes for 
a p-pulse controlled rectifier with a two-pole LC output 
filter (as in Chapter III) would be beneficial. 

5. A Better Solution for the Output Voltage Ripple 
for a Controlled Rectifier with a Two-Pole LC 
Output Filter 

In Chapter IV, approximations were introduced in the 
determination of the output voltage ripple for a controlled 
rectifier with a two-pole LC output filter. This had to do 

with the fact that |/ AC | was calculated from V AC and the 

overall load impedance |Z toto/ | only at the frequency pco . While 
po) is the dominant frequency, load impedances at frequency 
multiples of pco also exist. A more accurate estimate of 

|/ AC | can be obtained by using the principle of superposition 

to solve for \Iac\ using several values of |Z toto/ | at frequency 
multiples of pco and adding them together. While this will 
provide a better estimate of |/ AC | , which will in turn result 
in a more accurate value for V C(RMS) , the fact remains that 
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V C(RMS ) is not perfectly sinusoidal, and V c{pp _ ripple) and %V c{pp/avg) 

will still be approximations, although, somewhat better 
than before. It would still be interesting to determine the 
error caused by this approximation. 
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APPENDIX. MATLAB® CODE 


A. BASIC LINE-TO-LINE VOLTAGE WAVEFORM PLOT FOR A 
HYPOTHETICAL THREE-PHASE CIRCUIT 


clc 

clear 


f =60; 
w=2*pi*f; 

t=linspace(0, 3/pi, 1000); 

Vab = sind(w*t); 

Vbc = sind(w*t-120); 

Vca = sind(w*t-240); 

Vba=-Vab; 

Vcb=-Vbc; 

Vac=-Vca; 

hold on 

plot (w*t, Vab, 'b' , w*t, Vbc, 'g' , w*t, Vca, 'r' , 'LineWidth' , 2) 

axis([0,360,-1,1]) 

xlabel( 'angle (degrees)') 

ylabel( 'Voltage (Volts)') 

legend( 'V_a_b' , 'V_b_c' , 'V_c_a' ) 

grid on 

hold off 

B. LINE-TO-LINE VOLTAGE WAVEFORM PLOT SHOWING SIX PULSES 
FOR A HYPOTHETICAL THREE-PHASE CIRCUIT 


clc 

clear 

f=60 ; 
w=2*pi*f; 

t=linspace(0, 3/pi, 1000); 

Vab = sind(w*t); 

Vbc = sind(w*t-120); 

Vca = sind(w*t-240); 

Vba=-Vab; 

Vcb=-Vbc; 

Vac=-Vca; 

hold on 

plot (w*t, Vab, 'b' , w*t, Vbc, 'g' , w*t, Vca, 'r' , 'LineWidth', 2) 
plot (w*t, Vba, 'c' , w*t, Vcb, 'k' , w*t. Vac, 'm' , 'LineWidth', 2) 
axis([0,360,-1,1] ) 
xlabel( 'angle (degrees)') 
ylabel( 'Voltage (Volts)') 

legend( 'V_a_b' , 'V_b_c' , 'V_c_a' , 'V_b_a' , 'V_c_b' , 'V_a_c' ) 

grid on 
hold off 
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C. CRITICAL INDUCTANCE CALCULATION AND PLOT FOR A 
CONTROLLED 6-PULSE RECTIFIER 


clear 

clc 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)); % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
B = sind(lil); % constant to make math cleaner/neater 
alpha_u = 85; % upper limit in degrees of firing angle alpha 
num_points = 5e5; % number of points for alpha 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcritl = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))); % critical 
inductance 

% for angles greater than the boundary angle 
hold on 

plot(alpha_d, Lcritl, 'k-.', 'LineWidth' , 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 

theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % angle in degrees at 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
Lcrit2 = (R./(p*A*w.*cosd(alpha_d) )).*... 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))); 
plot(alpha_d, Lcrit2, 'k' , 'LineWidth', 2) 

% title('Critical Inductance vs. \alpha for a 6 Pulse Controlled 
Rectifier') 

xlabel( '\alpha (degrees)') 

ylabel( 'Critical Inductance (Henry )' ) 

legend (' \alpha \geq 10.08 / '\circ' , '/alpha \leq 10.08 / '\circ' , . . . 

'location', 'North') 
grid on 
hold off 
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D. CRITICAL INDUCTANCE CALCULATIONS AND PLOT FOR 
CONTROLLED 6-, 12-, 18-, 24-, AND 96-PULSE RECTIFIERS 


clear 

clc 

hold on 

alpha_u =30; % upper limit in degrees of firing angle alpha 
num_points = 5e5; % number of points for alpha 

%% 6 pulse 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)); % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
B = sind(lil); % constant to make math cleaner/neater 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcritl = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))) ; % critical 
inductance 

% for angles greater than the boundary angle 
plot(alpha_d, Lcritl, 'r-.', 'LineWidth' , 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 

theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % angle in degrees at 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
Lcrit2 = (R./(p*A*w.*cosd(alpha_d) 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))) ; 
plot(alpha_d, Lcrit2, 'r' , 'LineWidth', 2) 

%% 12 pulse 

p = 12; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)) ; % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
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B = sind(lil); % constant to make math cleaner/neater 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcrit3 = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))); % critical 
inductance 

% for angles greater than the boundary angle 
plot(alpha_d, Lcrit3, 'm -.' , 'LineWidth', 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 

theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % angle in degrees at 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
Lcrit4 = (R./(p*A*w.*cosd(alpha_d) )).*... 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))); 
plot(alpha_d, Lcrit4, 'm' , 'LineWidth', 2) 

%% 18 pulse 

p = 18; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)); % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
B = sind(lil); % constant to make math cleaner/neater 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcrit5 = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))); % critical 
inductance 

% for angles greater than the boundary angle 
plot(alpha_d, Lcrit5, 'g-.', 'LineWidth', 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
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angle in degrees at 


% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 
theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
Lcrit6 = (R./(p*A*w.*cosd(alpha_d) 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))); 
plot(alpha_d, Lcrit6, 'g' , 'LineWidth' , 2) 

%% 24 pulse 

p = 24; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)); % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
B = sind(lil); % constant to make math cleaner/neater 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcrit7 = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))); % critical 
inductance 

% for angles greater than the boundary angle 
plot(alpha_d, Lcrit7, 'k-.', 'LineWidth', 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 

theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % angle in degrees at 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
Lcrit8 = (R./(p*A*w.*cosd(alpha_d) )).*... 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))); 
plot(alpha_d, Lcrit8, 'k' , 'LineWidth', 2) 

%% 96 pulse 

p = 96; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
lil = 90-(360/(2*p)); % lil = lower integration limit for start of a 
pulse 

A = cosd(lil); % constant to make math cleaner/neater 
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B = sind(lil); % constant to make math cleaner/neater 

alpha_b_r = atan((p/pi)-(B/A)); % boundary angle alpha in radians that 

% defines transition point between submodes of continuous 
conduction 

alpha_b_d = alpha_b_r*180/pi; % boundary angle alpha in degrees 
w = 2*pi*60; % radian frequency for 60Hz 
R = 1; % unit resistance chosen for ease of comparison 

% code for angles greater than the boundary angle 

alpha_d = linspace(alpha_b_d, alpha_u, num_points); % range of alpha 
% in degrees for Lcritl 

Lcrit9 = ((R*tand(alpha_d))/w)*(1-((pi*B)/(p*A))); % critical 
inductance 

% for angles greater than the boundary angle 
plot(alpha_d, Lcrit9, 'b-.' , 'LineWidth', 2) 

% code for angles less than the boundary angle 

alpha_d = linspace(0, alpha_b_d, num_points); % range of alpha 
% in degrees for Lcrit2 
alpha_r = alpha_d*pi/180; % alpha in radians 

theta_d = asind((p*A*cosd(alpha_d))/pi) - lil; % angle in degrees at 
which 

% the inductor current will be minimum for this submode 
theta_r = theta_d*pi/180; % theta in radians 
LcritlO = (R./(p*A*w.*cosd(alpha_d))).* ... 

(pi.*cosd(theta_d+lil)+ ... 

(p*A.*sind(alpha_d))+ ... 

((p.*theta_r-p.*alpha_r-pi)*A.*cosd(alpha_d))); 
plot(alpha_d, LcritlO, 'b' , 'LineWidth', 2) 

%% Plot Formatting 


grid on 

%%title('Critical Inductance vs. \alpha for an n-Pulse Controlled 
Rectifier') 

xlabel( '\alpha (degrees)') 

ylabel( 'Critical Inductance (Henry )' ) 


legend( '6-pulse, 

\alpha 

\geq : 

< 

00 

o 

o 

\circ' , 

' 6-pulse, 

\alpha 

\leq 

10 

. 08 

A \circ' , ... 










'12-pulse, 

\alpha 

\geq 

5.01 

A \circ' , 

'12-pulse 

, \alph 

a \leq 

5. 

01 A 

\circ' , . . . 










'18-pulse, 

\alpha 

\geq 

3.34 

A \circ' , 

'18-pulse 

, \alph 

a \leq 

3. 

34 A 

\circ' , . . . 










'24-pulse, 

\alpha 

\geq 

2.50 A 

\circ' , 

'24-pulse, 

\alpha 

\leq 

2 . 

50 A 

\circ' , . . . 










'96-pulse, 

\alpha 

\geq 

0.63 A 

\circ' , 

'96-pulse, 

\alpha 

\leq 


0.63 A \circ' , . . . 

'location', 'Northwest') 
hold off 

E. VOLTAGE RIPPLE CALCULATIONS AND PLOT FOR A CONTROLLED 
6-PULSE RECTIFIER WITH A RESISTIVE LOAD 


clear 


100 


clc 


%% 6-pulse 


p = 6; % p = pulse count; 
Vm = 1; 

lil = 90-(180/p); % lil = 
uil = 90+(180/p); % uil = 
A = cosd(lil); % constant 
E = cosd(2*lil) ; 

F = sind(2*111); 


i.e., p = 6 for 6-pulse controlled rectifier 


lower integration limit for start of a pulse 
upper integration limit for start of a pulse 
to make math cleaner/neater 


alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 


% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl."2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
figure (1) 

plot(alpha_d, Vavgl, 'k' , . . . 
alpha_d, Vrmsl, 'r' ,... 
alpha_d, Vrms_ripplel, 'b' , ... 

alpha_d, Vpp_ripplel, 'm' , 'LineWidth' , 2) 
xlabel (' \alpha (degrees)') 
ylabel( 'Voltage (Volts)') 

legend( 'V_a_v_g' , 'V_R_M_S' , 'V_R_M_S_-_r_i_p_p_l_e' , 'V_p_p_- 

_r_i_p_p_l_e' , 'location', 'NorthEast') 

grid on 

hold on 

figure (2) 

plot(alpha_d, Vpp_by_Vavg_percentl, 'k' , 'LineWidth', 2) 
xlabel (' \alpha (degrees)') 
ylabel( '%V_p_p_/_a_v_g (%)' ) 
grid on 

axis([0 uil 0 1000]) 
hold on 


% code for boundary angle 1 < alpha < boundary angle 2 

alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 

alpha in degrees 

Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)); 
Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 
figure (1) 

plot(alpha_d, Vavg2, 'k' ,... 
alpha_d, Vrms2, 'r' ,... 
alpha_d, Vrms_ripple2, 'b' , ... 

alpha_d, Vpp_ripple2, 'm' , 'LineWidth', 2) 

101 


'k' , 


'LineWidth' , 2) 


figure (2) 

plot(alpha_d, Vpp_by_Vavg_percent2, 


range of alpha in 


% code for alpha > boundary angle 2 
alpha_d = linspace(alpha_b2_d, uil, num_points); 
degrees 

Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + lil)); 

Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+ ... 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/ (4*pi))); 
Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

Vpp_ripple3 = Vm*sind(alpha_d + 111); 

Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 
figure (1) 

Vavg3, 'k' , . . . 

Vrms3, 'r' , . . . 

Vrms_ripple3, 'b' , . . . 

Vpp_ripple3, 'm' , 'LineWidth', 2) 


plot(alpha_d, 
alpha_d, 
alpha_d, 
alpha_d, 
figure (2) 
plot(alpha_d. 


Vpp_by_Vavg_percent3, 'k' , 'LineWidth', 2) 


hold off 


F. VOLTAGE RIPPLE CALCULATIONS AND PLOT FOR A CONTROLLED 
6-, 12-, 18-, 24-, AND 96-PULSE RECTIFIER WITH A 

RESISTIVE LOAD 

clear 

clc 


hold on 
%% 6-pulse 


p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
Vm = 1; 


lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
E = cosd(2*lil) ; 

F = sind(2*111); 

alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 


% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl. A 2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
plot(alpha_d, Vpp_by_Vavg_percentl, 'r' , 'LineWidth', 2) 


% % code for boundary angle 1 < alpha < boundary angle 2 
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% alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 
alpha in degrees 

% Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

% Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 

% Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

% Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)) ; 

% Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent2, 'r—', 'LineWidth', 2) 

~o 

% % code for alpha > boundary angle 2 

% alpha_d = linspace(alpha_b2_d, uil, num_points); % range of alpha in 
degrees 

% Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + lil)); 

% Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+... 

"6 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/(4*pi))); 

% Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

% Vpp_ripple3 = Vm*sind(alpha_d + lil); 

% Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent3, 'r-.', 'LineWidth', 2) 

%% 12-pulse 

p = 12; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
Vm = 1; 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
E = cosd(2*111) ; 

F = sind(2*lil); 

alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 

% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl. A 2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
plot(alpha_d, Vpp_by_Vavg_percentl, 'm' , 'LineWidth', 2) 

% code for boundary angle 1 < alpha < boundary angle 2 

alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 

alpha in degrees 

Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)); 

Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 

plot(alpha_d, Vpp_by_Vavg_percent2, 'm— ' , 'LineWidth', 2) 

% % code for alpha > boundary angle 2 
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range of alpha in 


% alpha_d = linspace(alpha_b2_d, uil, num_points); % 
degrees 

% Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + 111)); 

% Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+... 

"6 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/(4*pi))); 

% Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

% Vpp_ripple3 = Vm*sind(alpha_d + lil); 

% Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent3, 'm-.', 'LineWidth', 2) 

%% 18-pulse 

p = 18; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
Vm = 1; 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
E = cosd(2*111) ; 

F = sind(2*lil); 

alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 

% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl. A 2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
plot(alpha_d, Vpp_by_Vavg_percentl, 'g' , 'LineWidth', 2) 

% code for boundary angle 1 < alpha < boundary angle 2 

alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 

alpha in degrees 

Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)); 

Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 

plot(alpha_d, Vpp_by_Vavg_percent2, 'g —' , 'LineWidth', 2) 

% % code for alpha > boundary angle 2 

% alpha_d = linspace(alpha_b2_d, uil, num_points); % range of alpha in 
degrees 

% Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + lil)); 

% Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+... 

'O 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/(4*pi))); 

% Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

% Vpp_ripple3 = Vm*sind(alpha_d + lil); 

% Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent3, 'g—.', 'LineWidth', 2) 
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%% 24-pulse 

p = 24; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
Vm = 1; 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
E = cosd(2*111) ; 

F = sind(2*lil); 

alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 

% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl. A 2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
plot(alpha_d, Vpp_by_Vavg_percentl, 'k' , 'LineWidth' , 2) 

% code for boundary angle 1 < alpha < boundary angle 2 

alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 

alpha in degrees 

Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)); 

Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 

plot(alpha_d, Vpp_by_Vavg_percent2, 'k— ', 'LineWidth', 2) 

% % code for alpha > boundary angle 2 

% alpha_d = linspace(alpha_b2_d, uil, num_points); % range of alpha in 
degrees 

% Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + lil)) ; 

% Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+... 

"6 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/(4*pi))); 

% Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

% Vpp_ripple3 = Vm*sind(alpha_d + lil); 

% Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent3, 'k—'LineWidth', 2) 

%% 96-pulse 

p = 96; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
Vm = 1; 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
E = cosd(2*lil) ; 
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F = sind(2*lil); 

alpha_bl_d = (uil-lil)/2; % first boundary angle in degrees 
alpha_b2_d = lil; % second boundary angle in degrees 
num_points = 1000; % number of points for alpha 

% code for alpha < boundary angle 1 

alpha_d = linspace(0, alpha_bl_d, num_points); % range of alpha in 
degrees 

Vavgl = (A*Vm*p*cosd(alpha_d))/pi; 

Vrmsl = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripplel = sqrt(Vrmsl. A 2 - Vavgl."2); 

Vpp_ripplel = Vm*(1-sind(alpha_d + uil)); 

Vpp_by_Vavg_percentl = (Vpp_ripplel ./ Vavgl)*100; 
plot(alpha_d, Vpp_by_Vavg_percentl, 'b' , 'LineWidth' , 2) 

% code for boundary angle 1 < alpha < boundary angle 2 

alpha_d = linspace(alpha_bl_d, alpha_b2_d, num_points); % range of 

alpha in degrees 

Vavg2 = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms2 = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple2 = sqrt(Vrms2. A 2 - Vavg2. A 2); 

Vpp_ripple2 = Vm*(sind(alpha_d + lil) - sind(alpha_d + uil)); 

Vpp_by_Vavg_percent2 = (Vpp_ripple2 ./ Vavg2)*100; 

plot(alpha_d, Vpp_by_Vavg_percent2, 'b —' , 'LineWidth', 2) 

% % code for alpha > boundary angle 2 

% alpha_d = linspace(alpha_b2_d, uil, num_points); % range of alpha in 
degrees 

% Vavg3 = ((Vm*p)/(2*pi))*(1 + cosd(alpha_d + lil)); 

% Vrms3 = (Vm/sqrt(2))*sqrt(((2+p)/4)-((alpha_d*p)/360)+... 

"6 

((E*p*sind(2*alpha_d))/(4*pi))+((F*p*cosd(2*alpha_d))/(4*pi))); 

% Vrms_ripple3 = sqrt(Vrms3. A 2 - Vavg3. A 2); 

% Vpp_ripple3 = Vm*sind(alpha_d + lil); 

% Vpp_by_Vavg_percent3 = (Vpp_ripple3 ./ Vavg3)*100; 

% plot(alpha_d, Vpp_by_Vavg_percent3, 'b—'LineWidth', 2) 

%% Plot Formatting 

grid on 

axis ( [0 30 0 60]) 

%title('%V_p_p_/_a_v_g % vs. \alpha for an n-Pulse Controlled Rectifier 

with a Purely Resistive Load') 

xlabel( '\alpha (degrees)') 

ylabel( '%V_p_p_/_a_v_g (%) ' ) 

legend (' 6-pulse, \alpha \leq 30 A \circ' , . . . 

'12-pulse, \alpha \leq 15 A \circ' , '12-pulse, 15 A \circ \leq 
\alpha \leq 75 A \circ' , . . . 

'18-pulse, \alpha \leq 7.5 A \circ', '18-pulse, 7.5 A \circ \leq 
\alpha \leq 80 A \circ' ,... 

'24-pulse, \alpha \leq 3.75 A \circ', '24-pulse, 3.75 A \circ \leq 
\alpha \leq 82.5 A \circ' , . . . 

'96-pulse, \alpha \leq 1.875 A \circ' , '96-pulse, 1.875 A \circ \leq 
\alpha \leq 88.125 A \circ' , . . . 

'location', 'Northwest') 
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% full legend below 

% legend('6-pulse, \alpha \leq 30 A \circ', '6-pulse, 30 A \circ \leq 
\alpha \leq 60 A \circ', '6-pulse, \alpha \geq 60 A \circ',... 

% '12-pulse, \alpha \leq 15 A \circ', '12-pulse, 15 A \circ \leq 

\alpha \leq 75 A \circ', '12-pulse, \alpha \geq 75 A \circ',... 

% '18-pulse, \alpha \leq 7.5 A \circ', '18-pulse, 7.5 A \circ \leq 

\alpha \leq 80 A \circ', '18-pulse, \alpha \geq 80 A \circ',... 

% '24-pulse, \alpha \leq 3.75 A \circ', '24-pulse, 3.75 A \circ \leq 

\alpha \leq 82.5 A \circ', '24-pulse, \alpha \geq 82.5 A \circ',... 

% '96-pulse, \alpha \leq 1.875 A \circ', '96-pulse, 1.875 A \circ 

\leq \alpha \leq 88.125 A \circ', '96-pulse, \alpha \geq 
88.125 A \circ', . . . 

% 'location', 'Northwest') 

hold off 

G. VOLTAGE RIPPLE CALCULATIONS AND PLOT FOR A CONTROLLED 
6-PULSE RECTIFIER WITH A TWO-POLE LC OUTPUT FILTER 


clear 

clc 

hold on 

%% initial constants 

alpha_u =30; % upper limit on firing angle alpha 
num_points = 10000; % number of points for alpha 
Vm = (100/sqrt(2))*sqrt(2); 

L0 = 25e-3; % inductance of load inductor 
CO = 50e-6; % capacitance of load capacitor 
R0 = 10; % resistance of load resistor 

LL = L0-(.5*L0); 

LH = L0+(.5*L0); 

CL = CO-(.5*C0); 

CH = C0+(.5*C0); 

%% 6p_L0_C0 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = L0; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 
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Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'k' , 'LineWidth' , 2) 

%% 6p_LL_C0 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = LL; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'r' , 'LineWidth', 2) 

%% 6p_LH_C0 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = LH; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 
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Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'b' , 'LineWidth', 2) 

%% 6p_L0_CL 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = L0; % inductance of load inductor 
C = CL; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'k —' , 'LineWidth', 2) 

%% 6p_L0_CH 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = L0; % inductance of load inductor 
C = CH; % capacitance of load capacitor 
R = R0; % resistance of load resistor 
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lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 

Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'm' , 'LineWidth' , 2) 

%% Plot Formatting 

grid on 

axis ( [0 30 0 20]) 

%title('Vpp-ripple/Vavg % vs. \alpha for a 6-Pulse Controlled Rectifier 

with a Two-Pole LC Output Filter') 

xlabel( '\alpha (degrees)') 

ylabel( '%V_p_p_/_V_a_v_g (%)' ) 

legend( 'L_0 C_0' ,... 

'L_L C_0' ,... 

'L_H C_0' ,... 

'L_0 C_L' ,... 

'L_0 C_H' ,... 

'location', 'Northwest') 
hold off 

H. VOLTAGE RIPPLE CALCULATIONS AND PLOT FOR A CONTROLLED 
6-, 12-, 18-, 24-, AND 96-PULSE RECTIFIER WITH A TWO- 

POLE LC OUTPUT FILTER 


clear 

clc 

hold on 

%% initial constants 

alpha_u =30; % upper limit on firing angle alpha 
num_points = 10000; % number of points for alpha 
Vm = (100/sqrt(2))*sqrt(2); 
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LO = 25e-3; % inductance of load inductor 
CO = 50e-6; % capacitance of load capacitor 
RO = 10; % resistance of load resistor 

%% 6p_L0_C0 

p = 6; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = LO; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = RO; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'r' , 'LineWidth' , 2) 

%% 12p_L0_C0 

p = 12; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = L0; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil); 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
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Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'm' , 'LineWidth' , 2) 

%% 18p_L0_C0 

p = 18; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = LO; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = RO; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'g' , 'LineWidth', 2) 

%% 24p_L0_C0 

p = 24; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = LO; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = RO; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil) ; 

Z_C = 1 / (j*120*p*pi*C); 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 
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alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))); 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'k' , 'LineWidth' , 2) 

%% 96p_L0_C0 

p = 96; % p = pulse count; i.e., p = 6 for 6-pulse controlled rectifier 
L = L0; % inductance of load inductor 
C = CO; % capacitance of load capacitor 
R = R0; % resistance of load resistor 

lil = 90-(180/p); % lil = lower integration limit for start of a pulse 
uil = 90+(180/p); % uil = upper integration limit for start of a pulse 
A = cosd(lil); % constant to make math cleaner/neater 
F = sind(2*lil); 

Z_C = 1 / (j*120*p*pi*C) ; 

Z_L = j*120*p*pi*L; 

Z_R_p_C = (Z_C * R)/(Z_C + R); 

Ztot = Z_L + Z_R_p_C; 
mag_Ztot = abs(Ztot); 

alpha_d = linspace(0, alpha_u, num_points); % range of alpha in degrees 
Vavg = (A*Vm*p*cosd(alpha_d))/pi; 

Vrms = (Vm/sqrt(2))*sqrt(1 + ((F*p*cosd(2*alpha_d))/(2*pi))) ; 
Vrms_ripple = sqrt(Vrms. A 2 - Vavg. A 2); 

Vac = Vrms_ripple; 
mag_Iac = Vac / mag_Ztot; 

VCrms = mag_Iac * abs(Z_C); 

Vpp_ripple = VCrms * 2 * sqrt(2); 

Vpp_by_Vavg_percent = (Vpp_ripple ./ Vavg)*100; 

plot(alpha_d, Vpp_by_Vavg_percent, 'b' , 'LineWidth', 2) 

%% Plot Formatting 


grid on 

axis ( [0 30 0 9]) 

%title('Vpp-ripple/Vavg % 
with a Two-Pole LC Output 
xlabel( '\alpha (degrees)') 
ylabel( '%V_p_p_/_V_a_v_g (%)') 
legend( '6-pulse' , ... 

'12-pulse' , . . . 
'18-pulse' , .. . 
'24-pulse' , . . . 


vs. \alpha for 
Filter') 


6-Pulse Controlled Rectifier 


r • • 
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'96-pulse' , 
'location' , 
hold off 


'Northwest' ) 
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